Search for Anomalous Production of Events with Two Photons and Additional 

Energetic Objects at CDF 

T. Aaltonen,^'^ J. Adelmanj^"^ B. Alvarez Gonzalez''/^ S. Amerio'''^,"''* D. Amidei,^^ A. Anastassov,^^ A. Annovi,^^ 
J. Antos,^^ G. Apollinari,^* A. Apresyan,''^ T. Arisawa,^* A. Artikov/'' J. Asaadi,^'^ W. Ashmanskas/^ 

A. Attal,-* A. Aurisano,^'* F. Azfar,'*^ W. Badgett/^ A. Barbaro-Galtieri,^^ V.E. Barnes,''^ B.A. Barnett,^^ 
P. Barria-''-'^,'^'^ P. Bartos/^ G. Bauer,33 P.-H. Beauchemin,34 F. Bedeschi,''^ D. Beecher,^! S. Behari,^^ 

G. Bellettini'^^'^^ J. Bellinger,6° D. Benjamin," A. Beretvas/^ A. Bhatti,^! M. Binkley,^* D. Bisello'*^,'''* 
I. Bizjak«,3i R E. Blair,^ C. Blocker/ B. Blunienfeld,^^ A. Bocci/'^ A. Bodek,50 V. Boisvert,^° D. Bortoletto,'^^ 

J. Boudreau,** A. Boveia/^ B. Brau°,^^ A. Bridgeman,^^ L. Brigliadori'^'^,^ C. Bromberg,^^ E. Brubaker,^^ 
J. Budagov,!^ H.S. Budd,^° S. Budd,^^ K. Burkett/^ G. Busetto'^'','''^ P. Bussey,^^ A. Buzatu,^^ K. L. Byruni,^ 

S. Cabrera^, 17 C. Calancha,^^ S. Camarda," M. Canipanelli,^^ M. Campbell,^^ F. Canellil^l« A. Canepa,46 

B. Carls,25 D. Carlsniitli,60 R. Carosi,^^ S. Carrillo",!^ S. Carron/^ B. Casal/^ M. Casarsa,!^ A. Castro^^^ 
P. Catastini/-'^,47 D. Cauz,^^ V. Cavaliere^/,'^^ M. Cavalli-Sforza," A. Cerri,^^ L. Cerrito',^! S.H. Cliang,^^ 

Y.C. Chen,! ^ Chertok,® G. Cliiarelli,^^ G. Chlacliidze,!^ F. Clilebana,!^ K. Cho,^^ D. Clioklieli,!^ J.P. Chou,^^ 
K. Chung°,i8 w.H. Chung,*^^ Y.S. Chung,^^ -p^ chwalek,^^ C.I. Ciobanu,^^ ^ a_ Ciocci/^,^^ A. Clark,^! D. Clark,^ 
G. Compostella,^* M.E. Convery,!^ J. Conway,^ M.Corbo,^^ M. Cordelli,^^ C.A. Cox,^ D.J. Cox,^ F. Crescioli'=^47 
C. Cuenca Almenar,^! J. Cuevas'',!^ R. Culbertson,!^ J.C. Cully,^^ D. Dagenhart,!** M. Datta,i** T. Davies,22 
P. de Barbaro,^° S. De Cecco,^^ ^ Deislier,^^ G. De Lorenzo,^ M. Del^Orso^^''7 C. Deluca,^ L. Demortier,^! 
J. Deng-^,17 M. Deninno,^ M. d'Errico'''^,^'' A. Di Canto^%47 G.P. di Giovanni,'*'^ B. Di Ruzza,^^ J.R. Dittniann,'^ 
M. D'Onofrio,* S. Donati'=^47 P. Dong,!^ T. Dorigo,'*^ S. Dube,^^ ^ Ebina,^** A. Elagin,^^ R. Erbacher,^ 
D. Errede,25 S. Errede,^^ N. Ersllaidat''^'*^ R. Eusebi,^^ H.C. Fang,^^ S. Farrington,43 W.T. Fedorko," R.G. Feild,^! 
M. Feindt,27 J.P. Fernandez,^^ C. FerrazzaSf,'*'^ R. Field,i^ G. Flanagan'', ''^ R. Forrest,^ M.J. Frank,^ M. Franklin,^^ 
J.C. Freenian,!^ I. Furic,!^ M. Gallinaro,^i J. Galyardt,!^ F. Garberson," J.E. Garcia,^! A.F. Garfinkel,''^ 
P. Garosi/^,'*^ H. Gerberich,^^ D. Gerdes,^^ A. Gessler,^^ S. Giagu'*'',^^ y. Giakoumopoulou,^ P. Giannetti,*^ 
K. Gibson,48 J.L. Ginimell,^^ CM. Ginsburg,!^ N. Giokaris,^ M. Giordani",^^ p. Giromini,20 M. Giunta,*^ 
G. Giurgm,26 V. Glagolev,!^ D. Glenzinski,i* M. Gold,^® N. Goldschmidt,!^ A. Golossanov,i* G. Gomez,!^ 
G. Gomez-Ceballos,33 M. Goncharov,^^ O. Gonzalez,^^ I. Gorelov,^*^ A.T. Goshaw,!^ K. Goulianos,^^ A. Gresele'''^,'"' 
S. Grinstein,^ C. Grosso-Pilcher,^"' R.C. Group, U. Grundler,^^ J. Guimaraes da Costa,^^ Z. Gunay-Unalan,^^ 

C. Haber,29 S.R. Hahn,i« E. Halkiadakis,^^ g _y. Han,50 J.Y. Han,^" F. Happacher,^^ K. Hara,^^ D. Hare,^^ 
M. Hare,^'^ R.F. Harr,^^ M. Hartz,^^ K. Hatakeyama,^ C. Hays,^^ jyj Heck,^^ J. Heinrich,46 M. Herndon,^" 

J. Heuser,27 S. Hewamanage,^ D. Hidas,^^ q q HilF," D. Hirschbuehl,^^ A. Hocker,!^ S. Hou,i M. Houlden,^^ 
S.-C. Hsu,29 R.E. Hughes,^" M. Hurwitz,!^ U. Husemann,^! M. Hussein,^^ J. Huston,^^ J. Incandela," G. Introzzi,'*'^ 
M. Iori'''',52 A. IvanovP,*^ E. James,!'^ D. Jang,!^ B. Jayatilaka," E.J. Jeon,^* M.K. Jha,^ S. Jindariani,i® 

W. Johnson,*^ M. Jones,''^ K.K. Joo,^^ S.Y. Jun," J.E. Jung,^^ T.R. Junk,i« T. Kamon,54 D. Kar,!^ 
P.E. Karchin,^9 Y. Kato™,42 R. Kepliart,i« W. Ketchum,!^ J. Keung,^^ V. Khotilovicli,54 B. Kilminster 
. Kim,28 H.S. Kim,28 H.W. Kim,^^ J.E. Kim,^^ M.J. Kim,^^ S.B. Kini,^^ S.H. Kim,^'^ Y.K. Kim," N. Kir 



, ^ N. Kimura,^® 

A.V. Kotwal,!^ M. Kreps,^^ 



tr.rj. rs^arcniu, i. r\.d,iu , rt. rtepuaru, vv . iveucnum, j. jxeung, v. rviioLuovien, d. i\ 
D.H. Kim,28 H.S. Kim,^^ H.W. Kim,^^ J.E. Kim,^^ M.J. Kim,^^ S.B. Kim,^^ S.H. Kim,^'^ Y.K. Kim, 

L. Kirsch,^ S. Klimenko,!^ K. Kondo,^® D.J. Kong,^^ J. Konigsberg,!^ A. Korytov,!^ A.V. Kotwal, 

J. Kroll,46 D. Krop,i4 N. Krumnack,^ M. Kruse,!^ V. Krutelyov," T. Kuhr,^^ N.P. Kulkarni,^^ M. Kurati 
S. Kwang," A.T. Laasanen,''^ S. Lami,^^ S. Lammel,^^ M. Lancaster, R.L. Lander,^ K. Lannon",^° A. La 
G. Latino/-'',^^ L Lazzizzera'''',^* T. LeCompte,^ E. Lee,^^ H.S. Lee," J.S. Lee,^^ S.W. Lee'",^^ S. Leone, 

.T n T,pwls 18 n _T T,i'n 29 T T.lnar-T-P 43 M T-lnHtrrpn 18 F, T,lnplps 46 A T,kfpr 21 F) O T.lfvlntspv 18 C. T.hi 



G. Latino/-'^,47 L Lazzizzera'''',^* T. LeCompte,2 E. Lee,54 H.S. Lee," J.S. Lee,28 S.W. Lee'",54 S. Leone,*^ 
J.D. Lewis,is C.-J. Lin,29 J. Linacre,^^ Lindgren,!^ E. Lipeles,^^ A. Lister,2i D.O. Litvintsev,!^ C. Liu,"*^ 
T. Liu,i^ N.S. Lockyer,46 A. Loginov,^! L. Lovas,!^ D. Lucchesi'''*,'*^ J. Lueck,27 p. Lujan,29 P. Lukens,!^ 
G. Lungu,^^ J. Lys,29 R. Lysak,i^ D. MacQueen,^'* R. Madrak,i^ K. Maeshima,!** K. Makhoul,^^ P. Maksimovic,26 

S. Malde,''^ S. Malik,3i G. Manca^^^ A. Manousakis-Katsikakis,^ F. Margaroli,^^ C. Marino, 2^ C.P. Marino, 2^^ 
A. Martin,^^ V. Martin'=,22 M. Martinez,^ R. Martmez-Ballarm,32 Mastrandrea,^^ M. Mathis,^^ M.E. Mattson,^^ 
P. Mazzanti,^ K.S. McFarland,^° P. Mclntyre,^^ R. McNulty^^o A. Mehta,^^ P. Melitala,24 A. Menzione,^^ 
C. Mesropian,^! T. Miao,i® D. Mietlicki,^^ N. Miladinovic,'^ R. Miller,^^ C. Mills,23 M. Milnik,27 A. Mitra,i 
G. Mitselmakher,i9 H. Miyake,^^ Moed,23 N. Moggi,^ M.N. Mondragon",!^ C.S. Moon,28 R. Moore,i*^ 
M.J. Morello,'*'^ J. Morlock,27 p. Movilla Fernandez, J. Miilmenstadt,29 A. Muklierjee,^® Th. Muller,27 P. Murat,^^ 

M. Mussini='=,6 J. Nachtman°,i® Y. Nagai,^^ J. Naganoma,^*^ K. Nakamura,^*^ I. Nakano,4i A. Napier,^'^ J. Nett,*^" 
C. Neu^4^ M.S. Neubauer,25 S. Neubauer,27 J. Nielsen9,29 L. Nodulman,2 M. Norman,io O. Norniella,25 E. Nurse,^! 
L. Oakes,43 g.H. Oh,!^ Y.D. 01i,28 L Oksuzian,!^ T. Okusawa,42 ^ Orava,24 K. Osterberg,24 S. Pagan Griso'*'^,^* 



2 



C. Pagliarone,^^ E. Palencia,^^ V. Papadimitriou,^^ A. Papaikonomou,^^ A. A. Paramanov,^ B. Parks,^° 

S. Pashapour,34 j_ Patrick,i« G. Pauletta",^^ M. Paulini,i=' C. Paus,33 T. Peiffcr,^^ D.E. Pcllctt,^ A. Pcnzo;"^-"^ 
T.J. Phillips," G. Piacentino,''^ E. Pianori ,"6 L. Pinera.^^ K. Pitts,^^ C. Plager,^ L. Pondroin,60 K. Potamianos,''^ 
O. Poukhov*/'' F. Prokoshinf/6 a. Pronko/^ F. Ptohos'/^ E. Pueschel/^ G. Punzi'^^ J. Pursley,^^ 
J. Rademacker'^,^3 A. Rahaman,^^ V. Ramakrishnan,^'^ N. Ranjan,^^ I. Redondo,^^ P. Renton,'*'^ M. Renz,^^ 
M. Rescigno,52 S. Ricliter,^^ F. Rimondi'^^'' L. Ristori ,"7 A. Robson,22 T. Rodrigo,!^ T. Rodriguez,^^ ^ Rogers,^^ 
S. Rolli,57 R. Roser,i8 M. Rossi,^^ R. Rossin," P. Roy,^^ A. Ruiz/^ j. Russ/^ V. Rusu/« B. Rutherford, 
H. Saarikko,24 A. Safonov,^^ W.K. Sakumoto,™ L. Saiiti",^^ L. Sartori,^^ K. Sato,-"^^ A. Savoy-Navarro,^^ 
P. Schlabach,!* A. Schniidt,^^ E.E. Schmidt,!^ M.A. Schmidt," M.P. Schmidt*, m. Schmitt,^^ T. Schwarz,^ 
L. Scodellaro,i2 A. Scribano-''^',^^ F. Scuri,*^ A. Sedov,^^ S. Seidel,^^ Y. Sciya,*^ A. Semenov,!^ L. Sexton-Kennedy, 
F. Sforza'=^4^ A. Sfyrla,^^ S.Z. Shalhout,'^^ T. Shears,^^ P.F. Shepard,^^ M. Shimojima*,^*^ S. Shiraishi," 
M. Shochct," Y. Shon,^" I. Shrcybcr,^'^ A. Simonenko,!^ P. Sincrvo,-^'' A. Sisakyan,!^ A.J. Slaughter,!^ 
J. Slaunwhite,40 K. Sliwa," J.R. Smith,^ F.D. Snider,!^ R. Snihur,^^ A. Soha,!^ S. Somalwar,53 V. Sorin,^ 
P. Squillacioti-'^/,*^ M. Stanitzki,*^^ R. St. Denis,^^ B. Stelzer,^^ O. Stelzer-Chilton,=^^ D. Stentz,^^ J. Strologas,^^ 
G.L. Stryckcr,3-'5 j g g^^j^ 28 a. Sukhanov,i« I. Siislov,!^ A. Taffard-'',^^ R. Takashima,'*! Y. Takciichi,''^^ r Tanaka,"! 
J. Tang," M. Tecchio,^^ P.K. Teng,i J. Thom'',!^ J. Thome," G.A. Thompson,^^ E. Thomson,''^ R Tipton,^! 
P. Ttito-Guzman,32 S. Tkaczyk.^^ D. Toback,^* S. Tokar,!^ K. Tollefson,^^ T. Tomura,^^ D. Tonelli,^^ 8. Torre,^^ 
D. Torretta," P. Totaro",'^'^ S. Tourncur,''^ M. Trovatoss,^^ S.-Y. Tsai,i Y. Tu,'*^ N. Turini-''/,^^ F. Ukegawa,^^ 
S. Uozumi,28 N. van Remo^tel^24 A. Varganov,^^ E. Vataga^s,-*^ F. Vazquez", G. Velev,!^ C. Vellidis,^ 
M. Vidal,32 I. Vila,i2 R. Vilar,!^ M. Vogel,^^ I. Volobouev'",^^ G. Volpi^^^^ P. Wagncr,^^ R.G. Wagner,^ 
R.L. Wagncr,i« W. Wagner'^'^,^^ J. Wagncr-Kuhr,^^ T. Wakisaka,42 R. Wallny,^ S.M. Wang,i A. Warburton,^^ 
D. Watcrs,3i j^j Weinberger,'^" J. Wcinclt,^^ W.C. Wester B. Whitchousc," D. Whitcson/,*^ A.B. Wicklund,^ 

E. Wicklund,!^ S. Wilbur," G. Williams,^" H.H. Williams,'"^ P. Wilson,i« B.L. Winer,*^ P. Wittich'\i« 
S. Wolbers,i8 C. Wolfe," H. Wolfe,^° T. Wright,^^ x. Wu,^! F. Wurthwein,!^ A. Yagil,i° K. Yamamoto,"^ 
J. Yamaoka,!^ U.K. Yang''," Y.C. Yang,^^ W.M. Yao,29 G.P. Ych,!^ K. Yi°,^^ J. Yoh,!^ K. Yorita,^^ T. Yoshida',42 
G.B. Yu,i7 I. Yu,28 S.S. Yu," J.C. Yun," A. Zanetti,^^ Y. Zeng,!^ X. Zhang,^^ Y. Zheng'',^ and S. Zucchelh'='=6 

(CDF Collaborationt) 

^Institute of Physics, Acadeinia Sinica, Taipei, Taiwan 11529, Republic of China 
^Argonne National Laboratory, Argonne, Illinois 60439 
^University of Athens, 157 71 Athens, Greece 
^Institut de Fisica d'Altes Energies, Universitat Autonoma de Barcelona, E-08193, Bellaterra (Barcelona), Spain 

^Baylor University, Waco, Texas 76798 
^Istituto Nazionale di Fisica Nucleare Bologna, '^'^ University of Bologna, 1-40127 Bologna, Italy 
^Brandeis University, Waltham, Massachusetts 02254 
^University of California, Davis, Davis, California 95616 
^University of California, Los Angeles, Los Angeles, California 90024 

University of California, San Diego, La Jolla, California 92093 
University of California, Santa Barbara, Santa Barbara, California 93106 
^^Instituto de Fisica de Cantabria, CSIC-University of Cantabria, 39005 Santander, Spain 
'^Carnegie Mellon University, Pittsburgh, PA 15213 
^"^Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637 
Comenius University, 842 48 Bratislava, Slovakia; Institute of Experimental Physics, O4O 01 Kosice, Slovakia 
Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 

^^Duke University, Durham, North Carolina 27708 
Fermi National Accelerator Laboratory, Batavia, Illinois 60510 
University of Florida, Gainesville, Florida 32611 
Laboratori Nazionali di Frascati, Istituto Nazionale di Fisica Nucleare, 1-00044 Frascati, Italy 
^''University of Geneva, CII-1211 Geneva 4, Switzerland 
^^Glasgow University, Glasgow G12 8QQ, United Kingdom 

^^Harvard University, Cambridge, Massachusetts 02138 
^^Division of High Energy Physics, Department of Physics, 
University of Helsinki and Helsinki Institute of Physics, FIN-OOOI4, Helsinki, Finland 
University of Illinois, Urbana, Illinois 61801 
^^The Johns Hopkins University, Baltimore, Maryland 21218 
^^Institut fiir Experimentelle Kemphysik, Karlsruhe Institute of Technology, D-76131 Karlsruhe, Germany 
'''^Center for High Energy Physics: Kyungpook National University, 
Daegu 702-701, Korea; Seoul National University, Seoul 151-742, 
Korea; Sungkyunkwan University, Suwon 440-746, 



3 



Korea; Korea Institute of Science and Technology Information, 
Daejeon 305-806, Korea; Chonnam National University, Gwangju 500-757, 
Korea; Chonbuk National University, Jeonju 561-756, Korea 
Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, California 94-720 
University of Liverpool, Liverpool L69 7ZE, United Kingdom 
University College London, London WCIE 6BT, United Kingdom 
^^Centro de Investigaciones Energeticas Medioambientales y Tecnologicas, E-28040 Madrid, Spain 
^^Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 

^'''Institute of Particle Physics: McGill University, Montreal, Quebec, 
Canada H3A 2T8; Simon Eraser University, Burnaby, British Columbia, 
Canada V5A 1S6; University of Toronto, Toronto, Ontario, 
Canada M5S 1A7; and TRIUMF, Vancouver, British Columbia, Canada V6T 2A3 
University of Michigan, Ann Arbor, Michigan 48109 
Michigan State University, East Lansing, Michigan 48824 
^^Institution for Theoretical and Experimental Physics, ITEP, Moscow 117259, Russia 
University of New Mexico, Albuquerque, New Mexico 87131 
^^Northwestern University, Evanston, Illinois 60208 
^"The Ohio State University, Columbus, Ohio 43210 
Okayama University, Okayama 700-8530, Japan 
Osaka City University, Osaka 588, Japan 
University of Oxford, Oxford 0X1 3RH, United Kingdom 
^"^Istituto Nazionale di Fisica Nucleare, Sezione di Padova-Trento, "^"^ University of Padova, 1-35131 Padova, Italy 
^^LPNHE, Universite Pierre et Mane Curie/IN2P3-CNRS, UMR7585, Pans, F-75252 Prance 
University of Pennsylvania, Philadelphia, Pennsylvania 19104 
^^Istituto Nazionale di Fisica Nucleare Pisa, University of Pisa, 
University of Siena and Scuola Normale Superiore, 1-56127 Pisa, Italy 
'^^ University of Pittsburgh, Pittsburgh, Pennsylvania 15260 
Purdue University, West Lafayette, Indiana 47907 
^'^ University of Rochester, Rochester, New York 14627 
^^The Rockefeller University, New York, New York 10021 
^^Istituto Nazionale di Fisica Nucleare, Sezione di Roma 1, 
Sapienza Universita di Roma, 1-00185 Roma, Italy 
^^Rutgers University, Piscataway, New Jersey 08855 
^'^ Texas A&M University, College Station, Texas 77843 
''^ Istituto Nazionale di Fisica Nucleare Trieste/Udine, 
I-34IOO Trieste, "University of Trieste/Udine, 1-33100 Udine, Italy 
^''University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
Tufts University, Medford, Massachusetts 02155 

^^Waseda University, Tokyo 169, Japan 
Wayne State University, Detroit, Michigan 48201 
University of Wisconsin, Madison, Wisconsin 53706 
Yale University, New Haven, Connecticut 06520 
(Dated: August 23, 2010) 

We present results of a search for anomalous production of two photons together with an electron, 
muon, T lepton, missing transverse energy, or jets using pp collision data from 1.1-2.0 fb^^ of 
integrated luminosity collected by the Collider Detector at Fermilab (CDF). The event yields and 
kinematic distributions are examined for signs of new physics without favoring a specific model of 
new physics. The results arc consistent with the standard model expectations. The search employs 
several new analysis techniques that significantly reduce instrumental backgrounds in channels with 
an electron and missing transverse energy. 
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I. INTRODUCTION 

Over the last twenty years, the rapid pace of develop- 
ments in phenomenology and model-building has left ex- 
perimentalists at the collider energy frontier with a wide 
array of new physics scenarios to investigate [l| . We are 
also assured that the number of models which have not 
yet been described is large. Since each search requires 
substantial resources, only a few new physics scenarios 
can be the focus of dedicated efforts. We address this 
problem by performing broad searches in available data 
samples for any discrepancy with the standard model 
(SM) [2[ in event yields or kinematic distributions. While 
this approach is not optimized for any particular scenario, 
it could possibly increase the chance of an unpredicted 
discovery. 

In this article we investigate a sample of data col- 
lected by the CDF II detector in pp collisions at ^/s = 
1.96 TeV at the Fermilab Tevatron. We restrict our- 
selves to a "baseline" sample with two isolated, cen- 
tral (0.05<|77|<1.05) photons (7) with £:t>13 GeV 8]. 
We then select subsamples which also contain at least 
one more energetic, isolated and well-identified object 
or where two photons are accompanied by large miss- 
ing transverse energy {fr). The additional object may 
be an electron (e), muon (/x), or r lepton (r). The 
is calculated from the imbalance in the energy of visible 
particles projected to the plane transverse to the beams. 
The integrated luminosity for each subsample varies from 
1.1 to 2.0 fb-i. 

The jj+X {X=e/^, r, and fJr) signatures are present 
in many new physics scenarios beyond the SM. Examples 
include models with the gauge-mediated supersymmetry 
breaking (GMSB) [|, extended Higgs sector [J|, techni- 
color models generation fermions ml , and theories 
with large extra dimensions 0]. 

The CDF collaboration has previously performed a 
search for anomalous production of two photons and an 
additional energetic object {fjT, e, /x, r, 7, jets, and b- 
quarks) in 85 pb~^ of the Tevatron Run I data @. Apart 
from the observation of a single eejj^ candidate event, 
the results were consistent with the SM predictions. The 
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eej^fpr event sparked considerable theoretical interest, 
because this signature is very rare in the SM and the 
event's topology is consistent with that of a decay of a 
pair of new heavy particles. In Run II, both the CDF [Toj 
and DO fll'] collaborations searched for production of 
77-1-^ events in the context of GMSB models using data 
corresponding to 0.20 fb"^ and 1.1 fb^^ of integrated lu- 
minosity, respectively. A search for anomalous produc- 
tion of 77-|-e//i events with energetic central photons and 
leptons {EY>2b GeV andp^>25 GeV/c) in 0.93 fb"! of 
data was performed by the CDF collaboration as a part 
of a broader signature-based search for new physics in 
I + ^ -\- X {I = e,fj, and X = jj, f}r) events Other 
signatures involving two photons were studied in CDF 
searches reported in Ref. [l3l - [T5| . The current model- 
independent analysis is improved upon previous diphoton 
searches both in terms of refined experimental techniques 
and of amount of data analyzed. It also probes a wider 
kinematic range compared to the CDF analyses reported 
in Ref. [12] and Ref. [15.]. 

This paper is organized as follows. It begins with a de- 
scription of the CDF II detector and the baseline dipho- 
ton sample. Then, each 77-l-X (X=e//z, r, and fr) 
subsample is discussed in separate sections where we de- 
scribe the definition of the subsamples, the calculation 
of the SM predictions, and the comparison of the data 
and the predictions. The details of several techniques 
are postponed to appendices. 



II. DETECTOR OVERVIEW 

The CDF II detector is a cylindrically symmetric ap- 
paratus designed to study pp collisions at the Fermilab 
Tevatron. The detector has been described in detail else- 
where [1^; only the detector components that are rel- 
evant to this analysis are briefly discussed below. The 
magnetic spectrometer consists of tracking devices inside 
a 3-m diameter, 5-m long superconducting solenoid mag- 
net which provides an axial magnetic field of 1.4 T. A set 
of silicon microstrip detectors (LOO, SVX, and ISL) [TtI - 
Hi] and a 3.1-m long drift chamber (COT) ^ with 96 
layers of sense wires measure momenta and trajectories 
(tracks) of charged particles in the pseudorapidity re- 
gions of \'q\<2 and 17?1<1 [1], respectively. Surrounding 
the magnet coil is the projective-tower-geometry sam- 
pling calorimeter, which is used to identify and measure 
the energy and position of photons, electrons, jets, and 
If^T- The calorimeter consists of lead-scintillator electro- 
magnetic and iron-scintillator hadron compartments and 
it is divided into a central barrel (]77]<1.1) and a pair 
of "end plugs" that cover the region 1.1<]77]<3.6. The 
central calorimeter is composed of towers with a segmen- 
tation of A77 X A(/) ~ 0.1 X 15°. The energy resolution 
of the central electromag netic calor imeter for electrons 
is a[ET)/ET = 13.5%/V£^T(GeV) © 1.5% [2l|, while 
the energy resolution of the central hadron calorimeter 
for charged pions that do not interact in the electromag- 
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nctic section is a{ET)/ET = 50% / y/B^^(GeV) ® 3% 
In the plug calorimeter, the segmentation varies from 
Ar] X A(j) ~ 0.1 X 7.5° for 1.1<|77|<1.8 to Atj x Acj) ~ 
0.6 X 15° for |77|=3.6. The corresponding plug elec- 
tromagnetic and hadron calorimeter energy resolutions 
are a{E)/E = 14.4%/ ^^^(GeV) © 0.7% and a{E)/E = 
74%/^£;(GeV) © 4%, respectively (H- The additional 
system in the central region is used for identification 
and precise position measurement of photons and elec- 
trons. Multiwire proportional chambers with cathode- 
strip readout (the CES system) are located at the depth 
of six radiation lengths (near shower maximum) in the 
central electromagnetic calorimeter. Cathode strips and 
anode wires, with a channel spacing between 1.5 cm and 
2 cm, running along the azimuthal (strips) and the beam 
line (wires) direction give location and two-dimensional 
profiles of the electromagnetic showers. The position res- 
olution of the CES is 2 mm for a 50 GeV photon. The 
electromagnetic compartments of the calorimeter are also 
used to measure the arrival time of particles deposit- 
ing energy in each tower (2^ . Muons from collisions as 
well as cosmic rays are identified using systems which 
are located outside the calorimeters: the central muon 
detector (CMU) and the central muon upgrade detec- 
tor (CMP) in the pseudorapidity region of |77|<0.6, and 
the central muon extension (CMX) for the pseudorapid- 
ity region of 0.6<|77|<1.0 ^51]. The CMU system uses 
four layers of planar drift chambers and detects muons 
with pr>1.4 GeV/c. The CMP system, located behind 
a 0.6 m thick steel absorber outside the magnetic re- 
turn yoke, consists of an additional four layers of planar 
drift chambers and detects muons with pT>2.2 GeV/c. 
The CMX detects muons with pT>lA GeV/c using four 
to eight layers of drift chambers, depending on the po- 
lar angle. A system of Cherenkov luminosity counters 
(CLC) f26j, located around the beam pipe and inside the 
plug calorimeters, is used to measure a number of in- 
elastic pp collisions per bunch crossing, and thereby the 
luminosity. 

The online event selection at CDF is done by a three- 
level trigger [275 system with each level providing a rate 
reduction sufficient to allow for processing at the next 
level with minimal deadtime. The Level- 1 uses custom 
designed hardware to find physics objects based on a sub- 
set of the detector information. The Level-2 trigger con- 
sists of custom hardware to do a limited event reconstruc- 
tion which can be processed in programmable processors. 
The Level-3 trigger uses the full detector information and 
consists of a farm of computers that reconstruct the data 
and apply selection criteria similar to the offline require- 
ments. 



III. DATA SELECTION AND EVENT 
RECONSTRUCTION 

The search for anomalous production of 77 + and 
77 -|- T events is performed with data corresponding to 



2.0±0.1 fb^^ of luminosity integrated from the beginning 
of Run II. The search for anomalous 77 + e/fi events 
utilizes a smaller dataset of l.lzbO.l fb~^ of integrated 
luminosity. The online (trigger) requirements and of- 
fline selection criteria that are common to all three final 
states are discussed below. Additional requirements of 
each analysis are explained separately. 

The inclusive 77 events are selected online by a three- 
level trigger that requires two isolated electromagnetic 
(EM) clusters with E^>12 GeV (diphoton-12 trigger) 
or two electromagnetic clusters with E^>18 GeV and 
no isolation requirement (diphoton-18 trigger). A de- 
tailed description of diphoton triggers can be found in 
Appendix lA II The triggered 77 candidate events are 
then subjected to the offline selection. Each event is re- 
quired to have two central EM clusters (photon candi- 
dates) inside a well-instrumented region of the calorime- 
ter (approximately 0.05<|ry|<1.05) with Et>13 GeV. 
For each photon candidate, the transverse shower pro- 
file in the CES and the amount of energy leaked into the 
hadron calorimeter must be consistent with those of a 
single electromagnetic shower. We distinguish photons 
from electrons by making sure that no high-p^ charged 
track points to the EM cluster. Both photon candidates 
are also required to be isolated energy clusters in the 
calorimeter in order to suppress background due to jets. 
More details of the standard photon identification criteria 
can be found in Appendix I A 21 To reduce contamination 
due to cosmic-ray, beam-related, and other non-collision 
backgrounds, the event must contain a well-reconstructed 
vertex, formed by tracks, with |z|<60 cm. If multiple ver- 
tices are reconstructed, the vertex with the largest J2 Pt 
of the associated tracks is selected. Unless noted other- 
wise, the transverse energy of all calorimeter objects is 
calculated with respect to this primary vertex. Finally, 
the arrival time of both photon candidates, corrected for 
average path length, has to be consistent with the pp 
collision time. It should be pointed out that due to the 
photon timing requirements, we are only sensitive to new 
physics processes where photons are produced in decays 
of new particles with small lifetime (<1 ns). 

Inclusive 77 events satisfying the above criteria form 
the baseline 77 signal sample used in all three analyses. 
This sample consists of real 77 even ts ( approximately 
30%), jet-j (45%) and jet-jet (25%) [ij events where 
one or both jets are misidentified as a photon. (An ob- 
ject misidentified as a photon is referred to as a "fake" 
photon.) The 77 -I- e//i, 77 -f- r, and 77 -I- frp candidate 
events are then selected from the base signal sample by 
requiring additional objects of interest or significant 
We also select a control sample of 77 events by applying 
less stringent photon identification requirements as dis- 
cussed in Appendix lA 21 To avoid an overlap with the 
signal sample, at least one photon candidate from the 
control sample must fail the standard photon cuts. The 
control sample is ideal for testing our analysis techniques 
because it has a similar event topology, but is dominated 
by background events (the fraction of real 77 events in 
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it is approximately 5%). 

Our baseline signal and control 77 samples consist 
of 31,116 and 42,708 events, respectively, in data cor- 
responding to 2.0 fb^^ of integrated luminosity. 



IV. SEARCHES FOR ANOMALOUS 
PRODUCTION OF 77 + X EVENTS 

In this Section, we describe in detail three separate 
searches for anomalous production of 77 + e//x, 77 + r, 
and 77 + ^ events. All analyses use the same baseline 
77 samples and utilize the same definitions of the addi- 
tional objects and kinematic variables: electrons, muons, 
T leptons, jets, soft unclustered energy, pr, and Ht- The 
Ht is defined as a scalar sum of IfjT and Et of all identi- 
fied photons, leptons, and jets. The detailed descriptions 
of these objects can be found in Appendices lA 3IIA 81 

A. The 77 + e/^ Final State 

We search for anomalous production of events contain- 
ing two photons and at least one additional electron or 
muon in data corresponding to l.liO.l fb"-"^ of integrated 
luminosity. The events of interest are derived from the 
77 baseline sample described in Section Hill The electron 
identification criteria are similar to those for the photon 
except that an electron candidate must have an energetic 
track pointing to the EM cluster. The momentum, p, of 
this track has to be consistent with the energy deposited 
in the EM calorimeter. The electron identification re- 
quirements are described in detail in Appendix IA3I 

A well-reconstructed COT track is identified as a muon 
candidate if it is matched to hit segments (stubs) in the 
central muon detectors, and its energy deposition pat- 
tern in the EM and HAD calorimeters is consistent with 
that left by a minimum ionizing particle. Details on the 
muon identification requirements can be found in Ap- 
pendix IA41 

The selected 776 and 77/.! events must have at least 
one electron or muon candidate with E^>20 GeV and 
p^>20 GeV/c, respectively. We compare the observed 
number of events and kinematic distributions in the data 
with those from our SM background predictions. Back- 
grounds for the 776 and 77/i signatures of new physics 
include: 

1. The SM production of Z— and W^lu in as- 
sociation with two photons (.^77, VF77), where 
photons are radiated from either the initial state 
quarks, charged electroweak boson (M^), or the 
final-state leptons. 

2. Backgrounds due to misidentified particles (fake 
photons or leptons) 

(a) electrons misidentified as photons (e.g., Z7 
events). 



(b) jets misidentified as photons (e.g, W^ + jet or 
Z7 -|- jet events), 

(c) jets misidentified as leptons (mostly 77 candi- 
date events with an additional jet). 

We describe below how these background contributions 
are estimated. 

The SM Z77 and M^77 contributions are estimated 
from Monte Carlo (MC) simulation. The MC samples are 
generated using the leading-order (LO) matrix-element 
generator madgraph [l^. The output of MADGRAPH 
is fed into pythia [l^l to carry out parton fragmenta- 
tion, simulation of the underlying event and additional 
pp interactions in the same bunch crossing, as well as 
initial- and final-state radiation. The output of pythia 
is then processed through the GEANT-based detector sim- 
ulation [sTl followed by the same reconstruction pro- 
gram as that for the data. To account for an imper- 
fect modeling of the CDF-II detector, the MC predic- 
tions are corrected for small differences (1-10%) in photon 
and lepton identification and trigger efficiences between 
data and MC. In addition, the LO cross sections pre- 
dicted by MADGRAPH are scaled to the next-to-leading- 
order (NLO) cross sections according to the i^-factors in 
Ref. [32]. These i^-factors are functions of the dilepton 
mass and Et of the highest-£'T photon, and their values 
range from 1.36 to 1.62 with an average of ^1.4 for the 
kinematic range of ^77 and 1^77 production. The un- 
certainty for this background prediction includes statisti- 
cal uncertainty due to the finite size of MC samples, 6% 
systematic uncertainty on the measured integrated lumi- 
nosity, and 7% systematic uncertainty on the .2^77 and 
M^77 cross sections due to uncertainties in the factoriza- 
tion and renormalization scales and parton distribution 
functions (PDF) [H. 

The rest of the background contains at least one 
misidentified object. Fake photons can arise from the 
hard bremsstrahlung of electrons in the detector mate- 
rial, inefficient electron-track reconstruction, or decays of 
tt'^, rp ^ or in jets. Although these background sources 
yield real photons in the final state, they are referred to 
as "fake" in this analysis, to be distinguished from the 
photons possibly produced by new exotic particles. The 
number of 77 + I {I = e, fi) events where at least one 
of the photons is faked by an electron is estimated from 
the l-f + e data (collected with the diphoton triggers). 
We obtain the prediction by applying a 6^7 misidenti- 
fication probability as a function of the Et of electron 
(about 2.7% and 1.5% for electrons with Et =20 GeV 
and 40 GeV, respectively) to the selected Ij + e events. 
More details about this misidentification probability and 
its uncertainty are included in Appendix IB 11 To esti- 
mate the number of 77 -|- Z events where at least one of 
the photons is a misidentified jet, we select the l"f + jet 
data collected with inclusive lepton triggers and multi- 
ply them by a jet^j misidentification probability as a 
function of the jet's Et (about 0.2% and 0.04% for jets 
with Et =13 GeV and >50 GeV, respectively). The de- 
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scription of the jet-^j misidentification probability and 
its associated uncertainty can be found in Appendix IB 21 
Also note that both + e and Z7 + jet samples may 
contain events with fake leptons. 

The last source of background is events with two real 
photons and a fake lepton from the direct diphoton pro- 
duction with additional jets. The number of "dipho- 
ton -|- fake lepton" events is obtained by applying £^7^- 
dependent misidentification probabilities from Ref. 
to the events with two photon candidates and an ob- 
ject which may fake a lepton. These objects are jets for 
electrons and isolated tracks for muons. The probabil- 
ity for a jet (isolated track) with Et{pt)—50 GeV to 
fake a central electron (muon) is ~0.01% (~1%). Details 
of the misidentification probabilities and their uncertain- 
ties are discussed in Appendix IB 31 According to earlier 
studies [2^ , only 29%±4% of observed diphotons are real 
diphoton events. In order to avoid duplication with the 
fake photon contribution estimated above, the number of 
"diphoton + fake lepton" events is multiplied by the real 
diphoton fraction (29%±4%), which gives the number of 
"real 77 -I- fake lepton" events. 

The fake photon signature can also be produced as a 
result of the bremsstrahlung of cosmic muons as they 
pass through the calorimeters. However, the probability 
for a real photon event to overlap with such a cosmic 
event is found to be very small: 1.5x 10~® (see Ref. [28|). 
Therefore, the cosmic backgrounds are negligible in the 
776 and 77^ searches. 

Table H] lists the expected and observed numbers of 
776 and 77/i events for > 13 GeV. At this stage 
of event selection, we observe three 776 events and zero 
77/i events. The leading background in the 776 chan- 
nel is due to events where at least one of the photons 
is a misidentified electron. The leading background in 
the 77/i channel is the electroweak production of .^77 
events. Figures [THU show several important kinematic 
distributions, including invariant mass, electron and pho- 
ton Et, pr, jet multiplicity, and Ht from data and the 
predicted backgrounds before applying the final selection, 
the silicon-track rejection (described next). 

TABLE I: Summary of the predicted and observed numbers 
of 77e and 77/t events before applying silicon-track rejection. 
The systematic uncertainty includes uncertainty due to MC 
statistics, uncertainties in the data luminosity, predicted cross 
sections, and the misidentification probabilities. 



Source 


electron 


muon 


^77 


0.90 


± 


0.09 


0.55 


± 


0.05 


Wj-y 


0.17 


± 


0.02 


0.09 


± 


0.01 


Ij + e— 7-7 


5.14 


± 


0.68 


0.02 


± 


0.02 


^7 + jet->7 


0.48 


± 


0.31 


0.13 


± 


0.09 


Fake Z-l-77 


0.13 


± 


0.05 


0.004 ± 0.004 


Total 


6.82 


± 


0.75 


0.79 


± 


0.11 


Data 


3 














TABLE II: Summary of the predicted and observed numbers 
of 776 and 77/1 events after applying silicon-track rejection. 
The systematic uncertainty includes uncertainty due to MC 
statistics, uncertainties in the data luminosity, predicted cross 
sections, and the misidentification probabilities. 



Source 


electron 


muon 


ZjJ 


0.82 


± 


0.08 


0.50 ± 0.05 


IF77 


0.15 


± 


0.02 


0.08 ± 0.01 


Z7 + e— >-7 


2.26 


± 


0.46 


0.004 ± 0.004 


Ij + jet^>7 


0.44 


± 


0.26 


0.12 ± 0.08 


Fake /+77 


0.12 


± 


0.05 


0.004 ± 0.004 


Total 


3.79 


± 


0.54 


0.71 ± 0.10 


Data 


1 










The dominant source of background in the 776 search 
are 7-1-66 events where one of the electrons is misrc- 
constructed as a photon. An electron may lose its track 
and be reconstructed as a photon because of catastrophic 
bremsstrahlung in the detector material in front of the 
COT. However, such an electron often leaves a few hits 
in the silicon detector and can be partially recovered 
by a special tracking algorithm (see Appendix lA 21 and 
Ref. [35.] for more details). We further compare the data 
and background prediction after removing events where 
at least one of the photons is matched to this type of 
electron track (silicon-track rejection). The silicon-track 
rejection suppresses ~80% of fake photons from the elec- 
tron bremsstrahlung for i?J>45 GeV (see Fig. [T^ from 
Appendix IB 31) while it has only ^1% inefficiency for real 
photons. Once this procedure is applied, the observed 
number of 776 events is reduced to one. The final back- 
ground predictions after the silicon-track rejection can be 
found in Table IE 

The robustness of our background estimation tech- 
nique is validated in the following three ways. First, we 
use an independent method to measure the misidenti- 
fication probabilities. Instead of counting electrons that 
satisfy standard central-electron criteria, we use the num- 
ber of electrons that satisfy photon-like electron criteria 
in the 6— >'7 fake rate denominator (see Appendix IB ip . 
The difference in the prediction of fake photons from 667 
events is ~4-ll%. Second, we cross-check if the 667 data 
(to which the 6— >7 fake rate is applied) contain signifi- 
cant number of fake electrons. We fit data with a com- 
bined likelihood of multiple electron identification vari- 
ables, where the signal shapes are obtained from electrons 
in Z decays and the background shapes are obtained from 
the sample enriched with fake electrons. The purity of 
electrons is estimated to be 97±2%. In addition, we com- 
pare the yields of 667 and iifij events in data and those 
predicted by Z7 madgraph MC. We divide the ratio of 
data to MC yields in the muon channel by the same ratio 
in the electron channel. If the 667 events contain signif- 
icant amount of fake electrons, the double ratio will be 
inconsistent with unity. The double ratio is found to be 
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FIG. 1: Kinematic distributions of the 776 events from the SM (dashed line) and total (solid line) background predictions as 
well as the three events observed in the data (marker). The total backgound includes SM and fake contributions. The gray 
boxes indicate the uncertainty in background determination. Each photon is required to have an _Et>13 GeV. Distributions 
from the top left to the bottom right are: a) three-body invariant mass; b) invariant mass of two photons; c) invariant mass 
of each electron-photon pair; d) Et of each photon; e) Et of the electron; f) Ht, scalar sum of fJr and Et of all identified 
photons, electrons, and jets; g) fJr; and h) number of jets with Et>15 GeV. 



1.10±0.15. Third, we examine the background estimate 
in larger samples, where either one photon and one elec- 
tron (67) or one photon, one electron, and one jet (ejj) 
events are required. The numbers of 67 and 67 j events 
in data are consistent with those from the background 
predictions within one standard deviation. 

To summarize, we do not observe any evidence for 
anomalous production of 776 and 77/i events. 



B. The 77 -f T Final State 



We search for events with two photons and a 
hadronically-decaying r lepton in data corresponding to 
2.0 fb~^ of integrated luminosity. These events are a 
subset of the baseline diphoton sample (see Section IIII[) 
with at least one r lepton candidate identified using the 
tight requirements and passing Et >15 GeV (see Ap- 
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FIG. 2: Kinematic distributions of the 776 events from the SM (dashed Hne) and total (sohd hne) background predictions. The 
total backgound includes SM and fake contributions. The gray boxes indicate the uncertainty in background determination. 
We observe zero events in the data. Each photon is required to have an Et>13 GeV. Distributions from the top left to the 
bottom right are: a) three-body invariant mass; b) invariant mass of two photons; c) invariant mass of each muon-photon pair; 
d) Et of each photon; e) pr of the muon; f) Ht, scalar sum of fJr and Et of all identified photons, muons, and jets; g) frr; 
and h) number of jets with _Et>15 GeV. 



pendix IXSj) . We select 34 77 + r candidate events. 

We consider two sources of backgrounds: the SM pro- 
duction of W—^Tv or Z-^TT with photons and 77 events 
with jets misidentified as t leptons. Other backgrounds 
are negligible. 

The electroweak backgrounds are estimated from Wj 
and MADGRAPH MC simulation. The LO order 
predictions are multiplied by the appropriate next-to- 
leading-order if-factors described in Section IIV Al and 



Ref. [11]. We find that these electroweak events with 
real leptons are dominated by events with at least one 
real photon, so we do not consider the case where both 
photons are misidentified jets. The simulation predicts 
the background from the cases of two real photons or one 
real with one fake photon to be 2.2±0.8 events, where the 
uncertainty comes from MC statistics. 

The dominant background in this search is from events 
with two reconstructed photons (which may be real or 
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FIG. 3: The mass of the two photons in 77 + r candidate 
events (marker) and the SM backgrounds (histogram). The 
gray boxes indicate the uncertainty in background determi- 
nation. 



fake), and jets, where one of the jets is misidentified as a 
r lepton. To estimate this background, we select events 
with two photons and a jet identified as "loose" r lep- 
ton candidate (see Appendix lA SP and apply the jet— >t 
misidentification probability (see Appendix IB 4[) . Since 
the misidentification probability is different for jets orig- 
inated by quarks or by gluons, and the ratio of quark jets 
to gluon jets here may be different from the one in the 
sample used to derive the jet—>-T misidentification prob- 
ability, we investigate a correction for the different types 
of jets in our sample. The probability for a quark jet to 
fake a r lepton is three times larger than the probabil- 
ity for a gluon jet. The process becomes more complex 
because a photon candidate may also be a misidentified 
jet, and the probability for a quark jet to fake a photon 
is ten times larger than for a gluon jet. We use pythia 
MC samples of diphotons, inclusive single photons and 
dijets to investigate the quark and gluon content of our 
data sample. Previous studies [1^ have determined that 
the baseline diphoton sample has approximately 30% real 
diphoton events, 45% events with a real photon and a jet 
misidentifed as a photon, and 25% events with two jets 
misidentified as photons. The simulations indicate that 
the quark-to-gluon ratio is significantly higher in the 
case of one real photon and one fake photon (80% quarks) 
than either of the other cases (approximately 30% in di- 
jet events and 40% in events with two real photons) and 
needs to be corrected for. We account for this effect by 
using two methods. In the first method, we simply apply 
the jet~^T misidentification rate and then make a cor- 
rection for the difference in the average quark-to-gluon 
ratio in the sample. In the second method, we allow for 
the possibility that quark jets will preferentially become 



misidentified photons, leaving the remaining jets to be- 
come misidentified taus. This method yields our reported 
central result, and the variation between the methods in- 
dicates a 13% systematic uncertainty which is added to 
the 20% systematic uncertainty in the misidentification 
probability. 



> 30 
o 

a 

in 25 



o 20 

UJ 

15 
10 
5 




10 



yy+T Search 
» Data, 2.0 fb"^ 
- Standard Model 



> 24F 
§22, 
in 20 - 
2 18^ 
« 16- 
i3 14- 
12^ 
10^ 
8 
6 

2 



20 30 40 50 60 70 80 90100 
Tau-lepton Ey (GeV) 



5 1015 20 25 30 35 40 45 50 



S! 35 
c 

I 30 
25 
20 
15 
10 
5 




5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 
Number of Jets, E^15 GeV 




50 100150 200250 300350400 
H,(GeV) 



FIG. 4: The kinematic distributions in 77-l-r candidate events 
(marker) and the SM backgrounds (histogram): a) Et of r 
lepton candidate; b) pr; c) number of jets with _Et>15 GeV; 
d) Ht, scalar sum of the transverse energies of photons, r 
lepton candidate, jets, and frp. The gray boxes indicate the 
uncertainty in background determination. 

The misidentified r background is 44±10 events and 
the total background estimate is 46±10 events, consis- 
tent with the 34 observed 77 -f r candidate events. We 
perform three checks of the methodology by predicting 
the size of 77 -I- r sample where the two photons are se- 
lected with the relaxed criteria (77 control sample de- 
scribed in Section IIIip or with one of the photons in 
the forward region (1.1<|?7|<2.0). The predictions for 
all control samples are consistent with the observations. 
Figures [3]|4] show several important kinematic distribu- 
tions for the selected 77 -f r candidate events and the 
predicted SM background. These distributions include 
the diphoton invariant mass, Et of a r lepton candidate, 
pr, jet multiplicity, and Ht- No excess is found above 
the SM background. 

In summary, we do not observe any evidence for the 
anomalous production of 77 -|- t events. 



C. The 77 -I- ^ Final State 

We search for the anomalous production of two pho- 
tons and large missing transverse energy (pr) h^ data 
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FIG. 5: The fh- distribution in all 77 candidate events from 
the baseline sample. The data (marker) is compared with the 
total background prediction (solid line with the gray band 
representing the total uncertainty). The total background 
prediction is a sum (shown by the stacked histograms) of 
the QCD, electroweak (dashed line), and non-collision (dash- 
dotted line) backgrounds. 



corresponding to 2.0 fb^^ of integrated luminosity. The 
subsample of 77 -I- frp events is derived from the base- 
line diphoton sample described previously in Section IlIII 
and in Appendix lA 21 The missing transverse energy is 
defined as an energy imbalance in the calorimeter (see 
detailed description of fir in Appendix lA 7p and it is 
an experimental signature of neutrinos or new particles 
that do not interact electromagnetically or strongly with 
the detector material. The fjTj however, can be mim- 
icked by a simple energy misreconstruction in SM events. 
Fluctuations in jet energy measurements are the most 
common source of such fake fr ■ Figure [5] shows the If^ 
distribution in the 77 baseline sample. This figure illus- 
trates that events with fake If^ are not only the domi- 
nant background in the region up to ^~40 GeV, but 
they also have a significant contribution even to the tail 
of pT distribution. A better separation between events 
with real and fake f^r can be achieved if a significance 
of the measured frp is considered rather than its abso- 
lute value. The ^-significance is a dimensionless quan- 
tity based on the energy resolution of jets and soft un- 
clustered particles. It also takes into account the event 
topology as shown in Appendix [Cl As it is demonstrated 
in Fig. [6l the ^-significance distributions have very dif- 
ferent shapes in events with fake and real frp'- exponen- 
tially falling (solid line) and almost fiat shapes, respec- 
tively. Thus, the ^-significance is an efficient tool in 
separating such events. We study 77 -I- pr events which 
pass three a priori ^-significance requirements: f}r- 
significance>3, 4, and 5. This choice of cut values has 



a straightforward motivation. If the 77 sample were only 
composed of events with fake fjr due to energy misrecon- 
struction in the calorimeter, then we would select 0.1%, 
0.01%, and 0.001% of the total number of events by re- 
quiring ^-significance>3, 4, and 5, respectively. On the 
other hand, studies with MC VF7— J-ei/ -|- 7 sample indi- 
cate that the ^-significance>3 (-^-significance>5) cut 
is ~100% (~90%) efficient for events with real |i^>35 
GeV (see Fig.©. 




^T-significance 

FIG. 6: The .^-significance distribution in all 77 candidate 
events from the baseline sample. The data (marker) is com- 
pared with the total background prediction (solid line with 
the gray band representing the total uncertainty). The to- 
tal background prediction is a sum (shown by the stacked 
histograms) of the QCD, electroweak (dashed line), and non- 
collision (dash-dotted line) backgrounds. The straight solid 
line represents the expected ^-significance distribution if all 
77 candidate events were to have fake frr due to the mea- 
surement fluctuations in the calorimeter (see Appendix [Cl for 
more details). 

We consider three major sources of background for 
the 77 + ^ signature: QCD (77, ^7, and jj where 
j^jet^jfake) events with large fake pr due to energy 
loss or mismeasurement in the calorimeter, electroweak 
(EWK) processes with real fir from neutrinos, and non- 
collision events with fake photons and Each of these 
sources is discussed below in the order of their impor- 
tance. All of the background estimation techniques are 
tested on a control sample of "loose" diphoton events 
described in Section IIIII and in Appendix lA 21 

Significant losses or fluctuations in energy measure- 
ments in the inclusive 77 sample, which can lead to 
considerable values of fake pr, happen only in a small 
fraction of events. However, the large production cross 
sections of QCD processes make them one of the largest 
backgrounds. We distinguish three types of QCD back- 
grounds: events with energy mismeasurement due to 
calorimeter energy resolution effects (QCD type-1); 77 



12 



> 

oio 

CO 

c 

0) 

> 

111 9 

10^ 



10 



10" 



— MC Wy-> ev+Y ^ 
■— after p-,.-significance>3.0 = 

after g-,.-significance>4.0- 

■ after p^-significance>5.ff| 




100 120 140 

?T(GeV) 



£T-significance>3.0 
gy-significance>4.0 
;y-significance>5.0" 




80 100 
(E^CGeV) 



FIG. 7: Figure a) demonstrates the fJr distributions in MC W^y-^eu+'y events before and after the .^-significance cuts. Figure 
b) illustrates the efficiency of the three ^-significance cuts for these events. 



candidate events with a wrong choice of the primary in- 
teraction vertex (QCD type-2); and 777 events where 
one of the photon candidates is lost in the calorimeter 
cracks (QCD type-3). 

The QCD type-1 background estimate is based on 
a ^-resolution model (metmodel) described in Ap- 
pendix O For each 77 data event, we generate ten 
pseudoexperiments to simulate fake fJr and calculate its 
significance given the event kinematics. In each pseudo- 
experiment, we smear the energies of jets and soft un- 
clustered particles using appropriate resolutions. The 
difference between the smeared and measured energy 
of the object is taken as its individual contribution to 
the total fake f/r- We predict the QCD background 
due to energy mismeasurements by counting the num- 
ber of pseudoexperiments that pass our ^-significance 
requirements divided by the number of pseudoexperi- 
ments per event. Kinematic distributions from these 
pseudoexperiments are then used as QCD background 
templates for data. The systematic uncertainty (23% 
for ^-significance>3, 47% for ^-significance>4, and 
130% for ^-significance>5) is evaluated by comparing 
the METMODEL expectations obtained with the default 
parameters to predictions obtained by varying each pa- 
rameter by one standard deviation (±1 cr). These param- 
eters and associated sources of systematic uncertainties 
are discussed in Appendix O The statistical and sys- 
tematic uncertainties are added in quadrature to obtain 
the total uncertainty. The predictions for the QCD type- 
1 background and their associated uncertainties can be 
found in Table HIl 

The background contribution due to 77 candidate 
events with mis-assigned primary vertex (QCD type-2) 
cannot be directly estimated by the metmodel because 



the energy resolution parameterization does not include 
this effect. The vertex misassignment occurs when a 77 
pair 36] is produced by a hard scattering interaction that 
overlaps with another interaction producing a vertex with 
the highest J^Pt of tracks. As a consequence of the 
wrong vertex choice, the Et of both photon candidates 
are incorrectly calculated, thus leading to fake frr- Al- 
though the effect is small, it can occasionally result in 
a large fake fr, for example, when two vertices are far 
apart and the photons are sufficiently energetic. We cor- 
rect for these mismeasurements by recalculating the Et 
of photons with respect to the vertex which gives the 
smallest value of fir- This procedure is verified to be 
valid for events with no intrinsic It is also tested in 
simulated 14^7— >'e^'7 events [3S] and data 67 events with 
fjT>2Q GeV. The selection of 67 events is discussed in 
Appendix lA 91 The effect is found to be small: after the 
procedure is applied, the number of simulated and data 
events with fJr > 20 GeV is reduced by 1% and 2%, re- 
spectively. In some fraction of events, however, the hard 
interaction completely fails to produce a reconstructed 
vertex and the vertex re-assignment cannot fix fake 
Since the metmodel cannot account for this contribu- 
tion, we employ a method based on a combination of 
data and Monte Carlo simulation to obtain the predic- 
tions. For this purpose, we use pythia 77 events [l^l 
passed through the detector simulation [3l[. These MC 
events also include additional interactions in the same 
bunch crossing that are modeled according to the lumi- 
nosity profile in data. We select only events where the 
hard scattering interaction resulting in a 77 pair does not 
produce a reconstructed vertex, and the primary vertex 
is created by tracks from an overlapping additional in- 
teraction. We will refer to such events as "no vertex" 77 
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events. The MC sample of "no vertex" events is normal- 
ized to the number of such events in real data (4.8±0.4% 
of the baseline 77 events). We then apply the standard 
analysis procedure to the sample and obtain the fraction 
of "no vertex" events in MC passing our ^-significance 
cuts. The systematic uncertainties on the QCD type-2 
background contribution include the MC statistical un- 
certainty (12%-24%), the uncertainty on the normaliza- 
tion factor (10%), the uncertainty due to the jet energy 
scale (7-8%), and the MC-data differences in the MET- 
MODEL parameterization (40%). The predictions for the 
QCD type-2 background and their associated uncertain- 
ties can be found in Table IIIII 

The 777 events are produced at a very low rate com- 
pared to that of 77 events. However a probability of 
losing a photon in calorimeter cracks is ^10%, so that 
the probability of losing one of the candidate photons 
in a tri-photon event is as large as ~30%. These events 
(QCD type-3) could reconstruct as 77 -|- ^ events. To 
reduce this background, we reject events if the ^ vector 
points along the direction (within |A(/)|<0.3) of a nar- 
row jet [39*1 located close to the calorimeter cracks at 
rj^O and It^I^I.I. The remaining contribution of the 
QCD type-3 events is estimated using a large inclusive 
PYTHIA 77 MC sample. We select reconstructed tri- 
photon events (i;J^'Sl3 GeV and E:^^>7 GeV) in MC 
and data. The numbers of reconstructed 777 candidates 
give us the MC-to-data normalization factor. To obtain 
an estimate of the remaining QCD type-3 background, 
we select PYTHIA tri-photon events at the generator level 
(before detector simulation) , apply the standard analysis 
procedure to these events, and multiply the result by the 
normalization factor described above. The systematic 
uncertainties for this background prediction is due to the 
following sources: 1) MC statistical uncertainty (24%- 
33%); 2) uncertainty on the normalization factor (19%); 
3) uncertainty due to MC-data differences in the met- 
MODEL parameterization (10%-44%); 4) jet energy scale 
uncertainty (10%-11%). The predictions for all sources 
of QCD backgrounds and their associated uncertainties 
can be found in Table IIIII 

Electroweak processes involving W^W and 
Z^vv It^t~ are the most common source of large 
real frr in collisions. There are three ways these 
processes can produce a 77 -f signature (listed in 
the order of importance) : 1 ) Wj and Zj events with 
one real and one fake photon; 2) VF77 and .^77 events 
where both photons are real; 3) W + jet, Z—>-t^t~, and 
Z + jet events where both photon candidates are fake 
photons. We estimate the EWK backgrounds by using 
W/Z + -f [38,] (for 1) and 2)) and inclusive W/Z (for 
3)) Monte Carlo events passed through the detector sim- 
ulation. We consider all three leptonic decay modes of 
W and Z bosons. To avoid an overlap between W/Z + 7 
and W/Z samples, we remove PYTHIA W/Z events where 
reconstructed photons are matched to generated photons 
originating from initial/final state radiation (ISR/FSR) 
of quarks or leptons. The MC-based predictions for 



the EWK backgrounds are then multiplied by a scale 
factor that diminishes possible data-MC differences 
and cancels out many of the systematic uncertainties 
(e.g., trigger efficiencies, acceptance and photon ID 
efficiencies, iiT-factors, modeling of ISR/FSR in MC, 
uncertainties in parton distribution functions, jet energy 
scale uncertainty, and luminosity uncertainty). This 
scale factor is obtained by comparing 67 -I- pr events 
(see Appendix lA 9P in data and MC. It is defined as the 
ratio of numbers of data and MC 67 events satisfying all 
analysis requirements. The resulting EWK background 
predictions and the corresponding uncertainties can be 
found in Table Hill The total uncertainties includes the 
MC statistical uncertainties (3.5-4.4%) and the MC- 
to-data normalization factor uncertainties (5.4-6.1%). 
The last uncertainty includes statistical uncertainties 
from data and MC 67 -I- fJr samples and systematic 
uncertainties associated with the purity of the e + 7 data 
sample and difference between the E/p-cnt efficiency 
(see Appendix IA9|) in data and MC. From Table |llTl 
one can see that the EWK processes are the dominant 
source of background when ^-significance>4. We find 
that 59-63% (30-40%) of the total EWK background 
for the 77 -|- pr signature comes from the electron 
(r lepton) decay channels of W and Z bosons. Note 
that the ^^-significance cuts are rather efficient for 
events with large real pr- for example, 84% and 68% 
of W+j^eiy+^ events pass the ^-significance>3 and 
_^-significance>5 requirements, respectively. 

The last remaining source of background is non- 
collision events where both photons and pr are fake. 
These events may either be caused by cosmic rays (CR) 
or beam halo (BH) muons depositing energy in the 
calorimeter. CR events are suppressed by requiring the 
EM timing of both photon candidates (Ti and T2) to 
be consistent with the collision time: |Ti,2|<6.7 ns and 
T2|<4.1 ns (more details are given in Appendix IA 2\i . 
BH events are removed by the topological cuts based on 
the distinct energy deposition pattern of BH muons trav- 
eling along the beam pipe. More details about CR and 
BH rejection cuts can be found in Appendix lA 21 The 
number of remaining BH events is estimated from the 
number of identified BH candidates and known rejection 
power of the BH cuts. The background contribution due 
to CR events is estimated based on the number of these 
events in the 30 ns< T12 <120 ns EM timing window 
and known efficiency of the cosmic rejection cuts (see 
Appendix lA 21 and Ref. [4l|). The total prediction for 
non-collision backgrounds can be found in Table IIIII The 
uncertainty for this estimate is dominated by the statis- 
tics in the samples of identified BH and CR events. 

The results of the search are presented in Ta- 
ble IIIII The total expected SM background for three 
^-significance cuts (-^-significance>3, 4, and 5) is 
71.7±7.5, 39.0±3.1, and 30.4±2.4 events, respectively. 
These predictions agree well with the observed numbers 
of data events: 82, 31, and 23. We also examine vari- 
ous kinematic distributions in data and SM backgrounds 
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TABLE III: The results of the search for anomalous production of 77 + fjr events. The data is compared to the background 
predictions for three values of the .^-significance cut. The quoted uncertainties include the effect of limited MC statistics as 
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/r-significance>3.0 ^ 


/r-significance>4.0 frr 


significance>5.0 


EWK 


35.4±2.2 


29.9±2.0 


25.9±1.9 


QCD type-1 
QCD type-2 
QCD type-3 


28.1±6.8 
4.4±2.0 
2.9±1.0 


3.6±1.8 
2.5±1.0 
2.2±1.0 


0.6±0.8 
1.5±0.7 
1.6±1.0 


Non-Collision 


0.9±0.3 
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0.8±0.3 


Total 


71.7±7.5 
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30.4±2.4 


Data 


82 


31 


23 



> 

C5 



c 

V 

> 



20 




18 




16 




14 




12 




10 




8 




6 




4 


r r 


2 




\ 




) : 



20 



YY+Ej Search a) 

• Data, 2fb"^ 
I I Uncertainty 
— QCD (fake 

EWK (real f^) 

ET-significance>3.0- 




40 



60 



80 100 

(GeV) 




50 100 150 200 
Leading Photon E,. (GeV) 




FIG. 8: The kinematic distributions for 77 + ^ candidate events with ^-significance>3: a) fr, the missing transverse energy; 
b) Et of leading photon candidate; c) Ht, the scalar sum of the transverse energies of photons, jets, and fr; and d) invariant 
mass, M-y^, of two photons . In all figures, the data (marker) is compared with the total background predictions (solid line with 
the gray band representing the total uncertainty). The total background prediction is a sum (shown by the stacked histograms) 
of the QCD and electroweak (dashed line) backgrounds. The non-collision background is too small to be visible on a plot with 
linear scale. 



for ^-significance>3 and ^-significance>5. Figures [5] distributions for the selected 77 -I- candidate events 
and [5] show the leading photon Et, Ht, and M^^ and the SM background predictions. Finally, Fig. [TUl il- 
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FIG. 9: The kinematic distributions for 77 + -^ candidate events with _^-significance>5: a) pr, the missing transverse energy; 
b) Et of leading photon candidate; c) Ht, the scalar sum of the transverse energies of photons, jets, and pr', and d) invariant 
mass, M^-y, of two photons. In all figures, the data (marker) is compared with the total background predictions (solid line with 
the gray band representing the total uncertainty). The total background prediction is a sum (shown by the stacked histograms) 
of the QCD and electroweak (dashed line) backgrounds. The non-collision background is too small to be visible on a plot with 
linear scale. 



lustrates multiplicities and Et distributions of extra jets 
and electrons in selected events. We observe good agree- 
ment between data and predicted background shapes for 
all studied kinematic distributions that are expected to 
be sensitive to production of new particles. 

In summary, we have searched for anomalous produc- 
tion of 77 -I- events in data corresponding to 2.0 fb^^ 
of integrated luminosity. No significant deviations from 
the SM background predictions are observed. 

V. CONCLUSIONS 

We performed a model-independent search for anoma- 
lous production of two photons with an electron, muon, 
T lepton, or large missing transverse energy. The anal- 



ysis of a jj+e/fi signature was performed using data 
corresponding to 1.1 fb~^ of integrated luminosity. Af- 
ter final selection, we observed one 77-l-e candidate event 
and zero 77+/^ events, in agreement with the expected 
background of 3.79±0.54 and 0.71±0.10 events, respec- 
tively. The kinematic properties of the 77-l-e event were 
consistent with the SM predictions. The silicon-track re- 
jection technique applied in this search allows for more 
than 60% reduction in the bremsstrahlung background 
(the dominant background in the electron channel) and 
has a promising potential for future searches with the 
j+e+X signature. 

The search for new physics in 77-l-r was based on data 
corresponding to 2.0 fb~^ of integrated luminosity. We 
observed 34 data events, in good agreement with the ex- 
pected background of 46±10 events. The kinematic dis- 
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FIG. 10: Properties of additional objects in 77 + f}r candidate events: a) the multiplicity of jets with _Et>15 GeV for events 
with _^-significance>3; b) the multiplicity of jets with _Et>15 GeV for events with _^-significance>5; c) the leading jet Et 
for events with _^-significance>3; and d) the electron Et for events with _^-significance>3. In all figures, the data (marker) 
is compared with the total background prediction (solid line with the gray band representing the total uncertainty) . The total 
background prediction is a sum (shown by the stacked histograms) of the QCD, electroweak (dashed line), and non-collision 
(dash-dotted line on figures a) and b) only) backgrounds. 



tributions of the selected events did not reveal any devi- 
ations from the SM predictions. 

The study of the ^"f+f}r signature was performed 
using data from 2.0 fb~^ of integrated luminosity. 
The events of interest were selected based on the 
significance, rather than a fixed -^-cut. This method 
proved to be very effective in rejecting events with fake 
pr, while remaining sensitive to new physics processes 
even with moderate values of pr (^^20-40 GeV). 
We selected 82, 31, and 23 data events with the pr- 
significance greater than 3, 4, and 5, respectively. These 
results are consistent with the expected SM background 
of 71.7±7.5, 39.0±3.1, and 30.4±2.4 events, respectively. 
The examined kinematic distributions for the observed 
events with ^-significance greater than three and five 



are in a good agreement with the predicted background 
shapes. The metmodel developed as part of the 77-1-^ 
search was also successfully applied to suppress multijet 
background with fake in the first observation of vec- 
tor boson pairs in a final state with two jets and pr at 



the Tevatron [42|. Finally, the reported in this paper 



model-independent analysis was later used as a basis for 
a search for supersymmetry with gauge-mediated break- 
ing in -f-f+pT events . The data samples used in these 
two analyses have a 60% overlap. 

In summary, no significant deviations from the stan- 
dard model were observed in the numbers of recorded 
events and their kinematic properties in signatures with 
two photons and an additional electron, muon, r lepton, 
or large frp. We also did not observe any new eejj-\-pr 
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candidate events, similar to the one reported in Ref. [9j. 
With improved analysis techniques and up to 20 times 
more data compared to the previous searches [ll-[ll|, this 
model-independent search is substantially more sensitive 
to new physics. 
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Appendix A: Definitions 
1. Diphoton Triggers 

There are two diphoton paths in the CDF three-level 
trigger: the first path requires two isolated electromag- 
netic (EM) clusters with £;t>12 GeV (diphoton-12), 
and the second path requires two electromagnetic clus- 
ters with Et>1S GeV and has no isolation requirement 
(diphoton- 18). The transverse energy of clusters is calcu- 
lated with respect to the nominal center of the detector 
at z=0 cm. The trigger requirements at each level are 
briefly described below. 

At Level-1, events with two towers with EM 
Et>8 GeV each are required. For each trigger tower, 
the amount of energy in the hadronic compartment of the 
calorimeter (E^^^) has to be consistent with that of an 
electromagnetic object. A trigger tower consists of two 
adjacent towers in the same calorimeter wedge, so that 
the granularity is approximately A77 x A(j> ~ 0.2 x 15°. 

The Level-2 requirements are different for the two trig- 
gers. The diphoton-12 trigger selects events if there are 
two isolated clusters (seeds) with EM Et>10 GeV each. 
The isolation (ISO) energy is calculated as a sum of the 
transverse energy in nine towers surrounding the seed 
tower according to five preset patterns. The ISO energy 
in each of the patterns has to be less than 3 GeV or 15% 
of the seed energy, whichever is larger. The diphoton-18 



trigger requires two towers with EM Et>16 GeV each 
at Level-2. 

The events are fully reconstructed at Level-3. At this 
level, for both triggers, the energy profile at the shower 
maximum (xcEs) '^^ each photon candidate has to be 
consistent with that of a single photon. The diphoton-12 
trigger selects events with two isolated photon candidates 
with i?T>12 GeV. The isolation energy at the level-3 is 
calculated as the sum of Et in all towers (except for pho- 
ton towers) within the cone of AR = y^Ajj^ + A<jP < 0.4 
centered around the photon candidate. This ISO energy 
has to be less than 2 GeV or 10% of the photon energy, 
whichever is larger. The diphoton-18 trigger has no iso- 
lation requirement and accepts events with two photon 
candidates with i<^T>18 GeV. Table [TVl gives a summary 
of all trigger requirements for events with EM objects in 
the central calorimeter. 



2. Photon Identification 

Photon candidates have to satisfy strict (also re- 
ferred to as "tight") photon identification requirements. 
The EM cluster has to be located inside the well- 
instrumented region of the CES chamber, away from the 
0-boundary of a calorimeter tower p^ . The energy de- 
position pattern in both transverse profiles at CES has 
to be consistent with that of an electromagnetic object. 
The ratio of the energy measured in the hadron (HAD) 
calorimeter to the EM energy, E^^^ /E^^, has to sat- 
isfy £;^^^/£;^*^<0.055-h0.00045x£:'^ requirement. To 
distinguish photons from electrons, no high-py charged 
track should point into the cluster {Ntrack<^ with track 
PT<l-0-l-0.005xii^T)- The main sources of " fake" photons 
are 7r° and rf produced in jets. These mesons are usu- 
ally produced in association with other particles. To re- 
duce this contamination from jets, the photon candidate 
must be isolated in the calorimeter and tracking cham- 
ber. To calculate the calorimeter isolation (cal-ISO), the 
Et deposited in the calorimeter towers within the cone of 
AR < 0.4 around the EM cluster is summed, and the Et 
due to the EM cluster is subtracted. The cal-ISO is then 
corrected for the photon's energy leakage into towers in 
the neighboring wedge and for the contribution from mul- 
tiple interactions in the same bunch crossing. The track 
isolation (track-ISO) is calculated as J^Pt of tracks in- 
side a cone Ai?<0.4 and satisfying \zyertex-ztrack\<5 cm. 
Both cal-ISO and track-ISO must be consistent with the 
amount of energy expected from the underlying event. 
In addition to calorimeter and tracking isolation, there 
should be no other significant energy {Et of 2"'' CES 
cluster) deposited in the CES chamber containing the 
photon candidate. Table |V] provides a summary of the 
photon identification requirements described above. 

We obtain the 77 control sample by selecting events 
where two photon candidates pass relaxed ( "loose" ) pho- 
ton identification requirements, but at least one of them 
fails the "tight" cuts. The main difference between 
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TABLE IV: Summary of the diphoton trigger requirements. 



Trigger Level 


Diplioton-12 


Diphoton-18 




EM Et> 8 GeV 


same 


Level- 1 


^HADi^EM ^ Q^^25 


same 




■^cluster 2 


same 




EM Et > 10 GeV 


EM Et > 16 GeV 


Level-2 


^HADi^EM ^ Q^25 


same 




E't?° < 3 GeV or E't^° /Et < 0.15 


N/A 




-^cluster 2 


same 




EM Et > 12 GeV 


EM Et > 18 GeV 


Level- 3 


^HAD i^EM ^ Q ^ 0.00045x_B/GeV if £ < 200 GeV 


same 




E't^° < 2 GeV or E^^°/Et < 0.1 


N/A 




shower profile; XCES < 20 


N/A 




Ncluster 2 


same 



TABLE V: Summary of the standard ("tight") and relaxed ("loose") photon identification requirements for the signal and 
control 77 samples, respectively. 



Cuts 


"Tight" photon ID 


"Loose" photon ID 




> 13 GeV 


same 


Shower profile in CES: 


<20 


same 




<0.055+0.00045x£;/GeV 


<0.125 


cal-ISO 


<0.1x_Bt if £'t<20 GeV or 


<0.15x£'t if £t<20 GeV or 




<2.0 GeV+0.02x(£T - 20 GeV) 


<3.0 GeV-F0.02x(£T - 20 GeV) 


track-ISO 


< 2.0 GeV+0.005 x Et 


< 5 GeV 


Ntracks in cluster 


< 1 


same 


track pt if Ntracks — 1 


<1.0 GeV+O.OOSx^T 


<0.25x£:t 


Et of 2"" CES 


<0.14x£;t if £t<18 GeV 


no cut 


cluster 


<2.4 GeV+O.Olx^T if £t>18 GeV 





"loose" and "tight" photon requirements is in the amount 
of allowed isolation energy (see Table |V]) . The resulting 
77 control sample is dominated by jet — 7 and jet — jet 
events where one or both photon candidates are faked by 
jets. The fraction of real 77 events in the control sample 
is only 5%. 

In addition to the photon identification requirements 
described above, we also apply cuts to remove contam- 
ination from non-collision sources. Muons produced in 
the beam halo are known to fake a photon signature [4l| . 
These energetic muons travel parallel to the beam pipe 
and deposit energy in many towers of one wedge, consis- 
tent with a minimum ionizing particle. When the muon 
undergoes energetic bremsstrahlung, it may also create 
one or two fake photon candidates. Probability for a sin- 
gle photon BH event to overlap with a collision event with 
a photon candidate is very low. Therefore, events with 
both fake photons from one muon are a dominant source 
of the BH background. We use this fact to suppress 
such BH events. We reject events if A(j)~^^ <0 .524 rad 
and if there are more than two hadronic and four cen- 
tral electromagnetic towers above O.l GeV threshold [isj . 
The efficiency of these requirements for collision events 
is estimated with data Z^e~^e~ events and found to be 
~100%. The rejection power of the cuts for beam halo 



events is found to be 90.4%±0.2%, as estimated using a 
very pure sample of beam halo events with two photon 
candidates located in the same calorimeter wedge. The 
criteria used to select this sample are discussed in detail 
in Ref. 

Muons from cosmic rays may also bremsstrahlung in 
the calorimeter and create fake photon candidates. To 
suppress these events, we use different approaches for 
data collected before and after the timing system in the 
EM calorimeter [2^j was installed. In the first 0.44 fb~^ 
of integrated luminosity for which the EM timing is not 
available, we reject events if there is a segment of hits 
("stub") in the muon drift chambers within a cone of 
30° around the direction of any of the photon candidates 
that is not linked to a track in the COT (trackless muon 
stub). This requirement rejects approximately 85% of 
cosmic rays and is approximately 98% efficient for 77 
events. In data from the later 1.6 fb^^ of integrated lu- 
minosity for which the EM timing is available, we reject 
events if one of the photon candidates has arrival time 
|T'i,2|>4crT or two photons have |Ar| = |Ti - r2|>4crAT, 
where (Tt=1-67 ns and (TAr=l-02 ns are the timing reso- 
lutions obtained by studying the EM timing of electrons 
from Z— >e+e~ events [46| . These EM timing require- 
ments reject 99.4% of cosmic rays, while they are 99.9% 
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efficient for prompt 77 events. 

Another source of fake pfiotons is electrons from elec- 
troweak processes wfiicfi are misreconstructed as prompt 
photons. This occurs when either an electron undergoes 
a catastrophic brcmsstrahhmg in the detector material 
in front of the COT or when its track does not get re- 
constructed. In both cases, electrons usually leave a few 
hits in the silicon detectors and their tracks can be par- 
tially recovered by a special tracking algorithm [ssj . This 
algorithm looks for silicon hits along two helix curves 
connecting vertex and EM cluster positions. The helix 
curvature is uniquely defined by the EM cluster Et, and 
two curves correspond to a positive and negative charge 
hypotheses. If any of the photon candidates is matched 
to such a track, we reject the event. This technique is 
used in the 77 -I- e/ fi and 77 -I- pr searches and it is re- 
ferred to as the "silicon track rejection" in the main text. 



3. Electron Identification 

We select electrons using the CDF standard criteria. 
An electron is characterized by a narrow shower in the 
electromagnetic calorimeter and a matching track (either 
in the COT or silicon detector) originating from the pri- 
mary vertex. The transverse EM energy, E^^, must be 
greater than 20 GeV. The £;HAD/^em ^^^^^^ ^las to be less 
than 0.055. The lateral energy distribution of the shower 
must be consistent with that for an electron. Candi- 
dates are required to be isolated in the calorimeter and 
to contain at least 90% of the total transverse energy 
within a cone of Ai?=0.4. For an electron detected in 
the central region (|77|<1.G5), the matching track is re- 
quired to be well reconstructed by the COT and have 
Pt>10 GeV/c. The ratio of the electron energy to track 
momentum, E/p, must be less than 2.0. Electrons from 
photon conversions are suppressed by rejecting the can- 
didates which have an oppositely-charged track with a 
small separation in the xy-plane and a minimal differ- 
ence in the polar angle. For an electron detected in the 
forward region (1.2<|?7|<2.0), the matching track is re- 
quired to have a minimum of three hits measured in the 
silicon detector. We do not apply further requirements on 
the forward matching tracks because fewer measurements 
per track are available and the momentum measured in 
the forward region is not as reliable as that measured in 
the central region. More details of electron identification 
can be found in Ref. [47| . 

4. Muon Identification 

We select muons using the CDF standard criteria. A 
muon is characterized by a well-reconstructed COT track 
which is matched to track segments (stubs) in the cen- 
tral muon detectors, and an energy deposition in the EM 
and HAD calorimeters consistent with a minimum ion- 
izing particle. The px, measured either with the COT 



only, or with the COT and the silicon detector if the sil- 
icon hits are available, must be greater than 20 GeV/c. 
Two types of muons are selected: CMUP (|?7|<0.6) and 
CMX (0.6<|r;|<1.0). The CMUP muon candidate re- 
quires a match between the track and the stubs in the 
CMU and CMP detectors. The CMX muon candidate 
requires a match to a muon stub in the CMX detector. 
In order to reduce the background from cosmic rays or 
hadrons which decay in fiight, we require the track to be 
consistent with originating along the beamline. Cosmic 
muons are further suppressed via their back-to-back track 
topology and asynchronous timing measured in the COT. 
More details on the muon identification can be found in 
Ref. 



5. Tau-Lepton Identification 

The r lepton has a ~18% branching fraction for decays 
into an electron or muon, with neutrinos. When this 
occurs, the event would be categorized in the e/fj, final 
state and addressed in that study. In the r lepton search 
we address only the hadronic decay modes. 

To identify the hadronic decays of r leptons , we re- 
quire a narrow cluster of one or three tracks and calorime- 
ter energy. This cluster must be consistent with a t lep- 
ton in several ways, inconsistent with an electron, and 
isolated from other nearby calorimeter energy. 

The clustering begins with a single tower with 
Et>Q GeV. Up to five more towers may be added to the 
cluster if they are adjacent and have i?T>l GeV. At least 
one high-quality track with pT>6 GeV/c must be associ- 
ated with the cluster. This track defines the origin point 
of the r lepton. The cone subtending an angle of Osig 
from the track direction defines the signal region where 
the T lepton decay products are expected. This angle 
is fixed to be 0.17 at low r lepton Et and is smaller for 
£:t>30 GeV, shrinking to 0.05 at Et^IOO GeV, allowing 
for greater rejection as the r lepton decay products be- 
come highly coUimated. A second cone given by 9<0.52 
defines an isolation annulus. The calorimeter and the 
shower maximum detector are used to define tt*' candi- 
dates in the r lepton signal cone and the isolation annu- 
lus. To reject electrons some hadronic energy, consistent 
with the observed signal-cone tracks, is required. 

The r lepton four-vector is defined by the total four- 
vector of the tracks and 7r° candidates in the signal cone. 
If the calorimeter cluster energy is significantly greater 
than this sum, the calorimeter cluster energy is used in- 
stead. The "visible" mass of the r lepton is found as the 
magnitude of this total four-vector. 

We define two levels of r lepton identification: a 
"loose" identification (used in studies and background 
techniques) and the standard or "tight" identification, 
used for the signal region search. 

Apart from the selection included in the reconstruction 
as described above, the loose identification requires only 
Et>15 GeV. The tight selection also requires visible t 
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lepton mass less than 1.8 GcV, total track px in the iso- 
lation cone less than 1.0 GeV, tt^ Et in the isolation cone 
less than 0.6 GeV, and one or three tracks in the signal 
cone, with total charge of ±1. 



6. Jets 

We reconstruct jets by using the cone clustering algo- 
rithm (42] with a cone radius R — -y/A^^ + A^ =0.4. 
Identified photons and electrons are removed from a list 
of jets. The jet energy is corrected for a non-linearity of 
the detector response and for contributions due to under- 
lying event and multiple interactions in the same bunch 
crossing [50|. Unless otherwise stated, we only consider 
jets with Et>1^ GeV and |r/|<3.0. 



9. The 67 Events in 77 + f}r Analysis 

We use inclusive 67 events to obtain a data/MC nor- 
malization factor for the MC-based estimate of the EWK 
backgrounds in the search for anomalous production of 
77 -|- IjtJT events. To minimize differences between 67 and 
77 samples, we obtain the 67 events in data and MC by 
using the same diphoton triggers (for data) and analy- 
sis selection procedures as used to derive our 77 baseline 
sample. In this selection, we treat an electron as a pho- 
ton (i.e., we apply the same cuts as in Table IV]) with only 
one exception that we also require the presence of a track 
pointing to an EM cluster. This track must satisfy the 
0.8<i?/p<1.2 requirement where E is the energy of the 
EM cluster and p is the track momentum. All additional 
tracks must pass the cuts listed in Table [V] 



7. Missing Transverse Energy and Ht 

The missing transverse energy, frr , is defined as an en- 
ergy imbalance in the calorimeter and it is a signature of 
neutrinos or new particles that do not interact with the 
detector material. The ^ is calculated from all calorime- 
ter towers with Et>0-1 GeV in the region |77|<3.6 accord- 
ing to fjT = — E'^fii where Hi is a unit vector that 
points from the interaction vertex to the i^^ calorimeter 
tower in the transverse plane. To improve resolution and 
reduce the number of events with large fake frp, we ap- 
ply corrections to the frp to account for a non-linearity of 
the detector response for jets with Et>15 GeV and for 
presence of reconstructed muons, which do not deposit 
their total energy in the calorimeter. 

One of the global kinematic characteristics of any hard 
scattering process is the total transverse energy of final 
products, Ht- We define Ht for each event as a sum of 
the transverse energies of all identified objects: photons, 
electrons, muons, visible energy of r leptons, jets, and 
pT- According to many theoretical models, new physics 
is expected to appear at large energy scales and may 
reveal itself in an anomalous rate of events with large 
values of Ht- 



8. Unclustered Energy 

The activity due to the underlying event and additional 
interactions in the same bunch crossing is characterized 

soft 

by the soft unclustered energy, Et - We calculate 

soft 

Et for each event by taking the difference between 
the total transverse energy in the event and transverse 
energies of all reconstructed photons, electrons, and jets: 

soft 



Appendix B: Fake Rates 
1. The e 7 Fake Rate 

Electrons may be misidentified as signal photons due 
to hard bremsstrahlung in the detector material, inef- 
ficiency of track reconstruction, or collinear final state 
radiation (FSR) . We measure the misidentification prob- 
ability, 'P{e—?''-f), using Drell-Yan Z/j* events. The 
7^(e— ^7) is defined as the ratio of the number of recon- 
structed Z/"f*^e^ events to the number of reconstructed 
Z/7*— T^ee events. The Et dependence of 7^(6-^7) is ob- 
tained from the simulation. The overall normalization 
is scaled by the ratio of data-to-MG probabilities mea- 
sured at £^T=40-50 GeV, around the Z peak. The elec- 
trons have been selected using two types of identifica- 
tion: 1) standard central-electron criteria in Section [A 31 
2) photon- like electron criteria in Section IXOl The result 
from identification- 1 has been used in the 776 and 77/i 
searches, while the result from identification- 2 has been 
used in the cross-checks of 776, 77/i, and jjfr searches. 
The V{e-^j) for identification- 1, before (7^(6—^7)^) and 
after (^^(e— >7)'^) applying silicon-track rejection, is mea- 
sured with the data and MG Drell-Yan samples and pa- 
rameterized as a function of the electron Et (in GeV) : 



-p^e ^ 7)^ = 5^ • (e-2.991-0.045.B^ 
^A oA 



0.007), (Bl) 

Vie ^ 7)^ = ■ Vie ^ 7)B • (1 - eiET)), (B2) 

where and are the data to MG scaling factors: 

= 1.08 ±0.09, 
= 1.32 ±0.20, 

and the eiET) is the efficiency of silicon-track rejection 
measured with the Drell-Yan MC: 



all 



eiET) = 0.405 • —= 



'dt 



(B3) 



'E% 



z = 0.086- (24.711 -^Jt) 
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FIG. 11: Probability for a CDF standard central electron to be misidentified as a standard central photon as measured in the 
Drell-Yan MC, before (a) and after (b) applying the silicon-track rejection. The misidentification probabilities (points) are 
parameterized as a function of electron Et at the parton level. The gray boxes indicate the systematic uncertainties in each 
Et bin due to the uncertainties on fit parameters. 
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FIG. 12: Efficiency for an electron reconstructed as a CDF 
standard central photon to be matched to the silicon-electron 
track, as a function of electron Et at the parton level, mea- 
sured in the Drell-Yan MC. The gray boxes indicate the sys- 
tematic uncertainties in each Et bin due to the uncertainties 
on fit parameters. 



The electron Et at the parton level is further trans- 
lated to the photon Et at the reconstruction level us- 



ing a simulation study. Figure [TT] shovifs the 7^(6— >7)^ 
and 7^(6— ?'7)^ measured in the Drell-Yan MC, without 
data-to-MC scaling factors applied. Figure [12] shows the 
^{Et)- The average 7'(e->7)^ is about 1.5% with an 11% 
fractional uncertainty, and the average 7^ (e—>-7)^ is about 
0.4% with a 17% fractional uncertainty. The uncertain- 
ties come from the limited size of Z data sample which 
determines the data-to-MC scaling factor, the variation 
of fitting methods which determines the number of Z 
candidates, and the difference between results measured 
in the diphoton and inclusive electron triggers. 



2. The jet 7 Fake Rate 

Hadrons in jets, such as 7r°, 77°, and K^, may decay 
into multiple photons. The segmentation of the CDF 
EM calorimeter is not sufficiently small to separate these 
photons and the standard reconstruction algorithm will 
reconstruct these hadron daughters as a single photon 
candidate. The probability to misidentify a jet as a signal 
photon, V{jet^^), has been measured in Ref. [s^, using 
data collected with inclusive jet triggers. The V{jet^^) 
is defined as the number of identified photon candidates 
times the fake-photon fraction (-Fqcd) and divided by the 
number of jets. The fraction -Fqcd is required because 
the identified photon candidates in the jet data will con- 
tain real photons not relevant to the fake rate. Ref. [s^l 
has determined Fqcd statistically by combining the fol- 
lowing information: (a) the lateral shower shape mea- 
sured in the wire and strip chamber (CES), (b) the extra 
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energy in a cone of Ai?=0.4 around the photon candidate 
(cal-ISO) measured in the calorimeters, and (c) the con- 
version rate measured in the central preshower detector. 
The V(jet—^^) is parameterized as a function of the Et 
of jet (in GeV) and found to be: 

V{]et ^ 7) = 10-3 . (e2.397-0.153.BT(jct) ^ q^q^^ (34) 

The fake photon Et is smaller than the original jet Et 
because the fake photon is often accompanied by other 
particles from that jet. The translation of the jet Et to 
the photon Et has been studied using simulations and 
is represented by a Gaussian distribution with a mean 
of 0.937 and a width of 0.048. The V{jet^-/) is about 
0.2% at =13 GeV and 0.04% for i;J>50 GeV, with a 
systematic uncertainty ranging from 50% to 200%. The 
sources of systematic uncertainties include: the differ- 
ences between the methods for determining i^qcD, the 
differences of results when using a loose photon candi- 
date as the fake denominator, variation of the mixture 
of quark jets and gluon jets, and variation of fragmenta- 
tion model in the simulation which changes the Gaussian 
function of Et translation. 

3. The jet e/fi Fake Rate 

Hadrons in jets may be misidentified as electrons due 
to inelastic charge exchange or the production of an ener- 
getic conversion electron. The inelastic charge exchange 
in the EM calorimeter, 7r~p—>-Tr^n or Tr+n— >-7r"p, results 
in a track in the COT due to the tt^ and an EM shower 
in the calorimeter due to the photons from tt^ decay. 
The combination of a charged track and an EM shower 
gives a fake electron candidate. Hadrons can also de- 
cay into muons before interacting with calorimeter (e.g., 
if"*"— j'^+f^) or pass through the calorimeter into the 
muon chamber (punch-through) with minimal interac- 
tion and give fake muon candidates. The probability to 
misidentify a jet as an electron or a muon, V{jet^e, fi), 
has been measured in Ref. [131, using data collected 
with inclusive jet triggers, ^^{jet^e, fi) is defined as 
the ratio of the number of identified electron/muon can- 
didates to the number of "fakeable" objects (denomi- 
nator). The fakeable object is a jet with uncorrected 
Et>4: GeV for central electrons, a jet with uncorrected 
Et>15 GeV for forward electrons, and an isolated track 
with pt>4 GeV/c and minimal extra energy in the cone 
of Ai?=0.4 for muons. A track is considered to be iso- 
lated if the total Et of calorimeter towers within the cone 
of Ai?<0.4 around the track is less than 4 GeV or less 
than 10% of track's momentum. The misidentification 
probabilities, parameterized as a function of the jet Et 
(in GeV), for the electrons are: 

Vijet ^ e™""'^') = 0.00013 + e-7-940-o.i94.ET(jct)^ ^g^^ 

Vijet 6^°™^''^) = 0.00032 + 0.000012 • £:T(jet).(B6) 

The translation of the jet Et to the electron Et has 
been studied using simulations. The ratio of electron Et 




Tau-lepton (GeV) 

FIG. 13: The solid line represents the probability for objects 
passing loose r lepton ID cuts to also pass tight r lepton 
ID cuts {jet — >■ r fake rate) as a function of r lepton Et, 
with the overlapping regions removed. The dashed lines are 
the systematic uncertainties on the jet — > r fake rate (±1 
standard deviation). 

to jet Et is represented by a Gaussian distribution with 
a mean of 0.89 and a width of 0.06. The misidentification 
probabilities, parameterized as a function of the track pT 
(in GeV/c), for the muons are: 

P{track = 0.00086 -I- 0.00017 • pTi^i), (B7) 

Vitrack /i^^^) = 0.00082 + 0.00020 • pt(m)- (B8) 

Since the misidentification probabilities were measured 
up to Et=50 GeV and pt=50 GeV/c and the misidenti- 
fication probabilities are expected to reach a plateau, we 
assign a constant value to all misidentification probabili- 
ties for £;t>50 GeV andpT>50 GeV/c. The V{jet^e, ji) 
averages saO.01% for central electrons, saO.04% for for- 
ward electrons, and «1.0% for central muons, with a 50% 
systematic uncertainty estimate provided by the varia- 
tion of results measured in different jet triggers. 



4. The jet t Fake Rate 

The probability of a quark or gluon jet to be misrecon- 
structed as an hadronically-decaying r lepton is mea- 
sured and then applied to a sample of jets to estimate 
the number of fake r leptons we expect in that sample. 

We measure this misidentification rate in a sample of 
inclusive jet triggers [sij , using only the energy clusters 
for which the trigger is fully efficient. This jet sample 
has a negligible fraction of real r leptons because the 
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rates for W/Z^t+X and c/b^r+X processes are very 
small compared to the jet production rates. Therefore, 
the measurement of the misidentification rate is straight- 
forward. We identify all loose and tight r leptons (See 
Section lASp and measure the rate by dividing the num- 
ber of tight T leptons by the number of loose r leptons 
as a function of the r-candidate Et- 

We check the misidentification rate by using it to pre- 
dict various distributions in the jet samples. We compare 
the number of r leptons observed and predicted as a func- 
tion of the first, second and third jet Et, the event total 
energy, the underlying event energy, the number of jets, 
the number of interactions in the event, and the distance 
the nearest jet. The only notable discrepancy is the case 
where the r lepton is close to a second jet, where the jet's 
energy tends to spoil the r lepton's isolation and reduce 
the fake rate. We include this effect in the application 
of the misidentification rate. The assigned systematic 
uncertainty of 20% accounts for any other discrepancies. 
The resulting function is shown in Fig. 1131 

Finally, we identify the primary source of each jet in 
inclusive PYTHIA MC jet samples by searching for the 
highest-E'r parton consistent with the jet direction. We 
then measure the misidentification rate in quark and 
gluon jets separately. We find the rate for gluon jets 
is approximately three times smaller than the rate for 
quark jets, which have a higher probability to fragment 
to a few energetic particles. 



Appendix C: The fJr Resolution Model 

A major sources of background in the 77 -I- ^ final 
state is diphoton candidate events with significant fake 
frr due to energy mismeasurement in the calorimeter. 
Given the large production rates for QCD processes (77, 
7-jei, and jet- jet), fluctuations in energy measurements 
can result in a considerable fraction of such events. We 
predict the shape of this fake fir and calculate its signif- 
icance on an event-by-event basis by means of the frp 
resolution model denoted as metmodel. 

The METMODEL is based on a simple assumption that 
fluctuations in energy measurements of jets, soft unclus- 
tered particles from the underlying event, and multiple 
interactions are the dominant sources of fake frr ■ There- 
fore, the individual contributions of each of these compo- 
nents to fake frr can be modeled, on average, by smearing 
their energies according to the corresponding energy reso- 
lution functions. Jets are the dominant source of fake frp 
because they are collimated sprays of energetic particles 
in a certain direction and may have large measurement 
fluctuations in that direction. The unclustered energy, 
on the other hand, tends to be uniformly spread in the 
calorimeter. Therefore, the portion of fJr due to this 
source is usually small and mostly results in a smear- 
ing of the jet component of fake Taking into account 
the above considerations and for reasons of simplicity, we 
model only the fake fJr due to mismeasurements of jets 
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FIG. 14: Example of the fake fir parameterization due to 
unclustered energy. Figure a) shows a two-Gaussian fit of the 
f}r^ distribution for pythia 77 events from one of the bins 
in \/^Et- Figure b) demonstrates how a width, a, of the 
leading Gaussian depends on the ^J^^Et- On both plots, 
points are pythia data and curves are the fit functions. 



and all soft unclustered energy (rather than individual 
unclustered particles) . 

The fjT resolution due to the soft unclustered energy is 
studied in the 77 control sample (see Appendix IA2p and 
Z/7*^e+e- events with 85 GeV/c2<Mee<97.5 GeV/c^. 
We fit distributions of x and y components of the for 
events without jets, Njet{ET>l5 GeV)=0, in small bins 
of -^/^ Et, with a sum of two Gaussian distributions. 
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We assume that both Gaussian distributions have the 
same mean, but different widths (cr and scale xct, respee- 
tively). From the individual fits of pr^ and fr^ distribu- 
tions, we obtain the mean, a, scale, and relative normal- 
ization, Norm, of two Gaussians for each bin of Et 
(bin size is 2 GeVs). The parameters are then fitted by 
simple polynomial functions of z—^J^ Et'- Po + Piz for 
a and the scale, po +Piz'^ for the mean, and Norm — po- 
These functions provide a parameterization of the unclus- 
tered energy contribution into the x and y components 
of the fake ^ in the event. The default set of param- 
eters is obtained from the 77 control sample. We also 
use the results of fits in the data Z— 7>e+e~ sample as 
an alternative set of parameters to study the associated 
systematic uncertainties. Figure [T4l demonstrates an ex- 
ample of the fJr^ resolution parameterization due to the 
unclustered energy in PYTHIA 77 events. Distributions 
for both X and y components of fJr look essentially iden- 
tical to those shown in Fig. [M) We also do not observe 
any large difference in the parameterization of the ^ res- 
olution due to unclustered energy between Z/j*~^e'^e~ 
and "loose" 77 events in data as well as between data 
and MC. 

To account for contributions from jets into the fake 
we obtain the jet energy resolution, JER, as a function of 
jet energy and pseudorapidity, E and 77. For this purpose, 
we use PYTHIA samples of dijet and jet events passed 
through the GEANT-based detector simulation. In these 
events, we reconstruct jets before {hadron jet) and after 
{detector let) the detector simulation by using the same 
cone clustering algorithm at both levels. The jet energy 
resolution is then defined as a ratio of the detector [E'^'^^) 
and hadron level {E'^°-'^) jet energies, J ER=E'^''^ / E'""^-l, 
for hadron and detector jets with pt>3 GeV/c that are 
matched within a cone of r])<G.l. Unlike the energy 
balance in dijet and Z-jet events, this definition of JER 
is mostly sensitive to detector effects and allows us to sig- 
nificantly minimize the dependence of resolution on the 
effects of initial and final state radiation. However, we 
still compare the dijet and Z-jet balance in data and MC 
to make sure that the simulation adequately describes the 
resolution. We fill JER histograms for jets in 5 GeV bins 
in jet energy and A77=0.2 bins in pseudorapidity. We fit 
these histograms by a linear combination of Gaussian and 
Landau functions of x, where x=—JER/{l + JER) en- 
sures stable fits in the entire range of jet energies. Exam- 
ples of fits for one particular 77-bin can be found in Fig. 1151 
These plots illustrate that our fit function successfully 
describes the jet energy resolution in a wide range of jet 
energies. It is also important to mention that the same 
fit function is used for all yy-bins. From the individual fits 
for each (i?-''^*,77)-bin, we obtain a relative normalization, 
C, and parameters of a Gaussian (mean and a) and Lan- 
dau (mean and a) fits. These parameters are plotted as 
a function of E^^* f or each ?7-b in, and fit with the follow- 
ing functions: a—y/po/E + pi; mean=po + piE + P2/ E; 

and C={po + pi^/E)/E + p2- This provides a smooth 
parameterization of JER for all reconstructed jets with 



i;j'=*>3 GeV and |77|<3.6. 
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FIG. 15: Examples of jet energy resolution (JER) fits using 
a linear combination of Gaussian and Landau functions of 
x^-JER/{l + JER) where JER = e'^"'' / E'^"'^ - 1 for two 
different jet energy bins: a) 20 GeV<E''^^<25 GeV and b) 
400 GeV<£"''='<405 GeV. 

We predict the shape of fake f!r based on the energy 
resolution functions described above. For each event, we 
produce a probability distribution function, V{fir), of all 
possible values of the fake fJr by smearing the energies 
of jets and unclustered energy according to these objects 
individual resolution functions in a large number of pseu- 
doexperiments. Then, we sum up these individual V{frr) 
distributions for all events to obtain a shape of the pre- 
dicted fake pT due energy mismeasurcmcnts in our data 



25 



sample. Technical details of how we generate V{frr) are 
given below. An example of this Vifr) distribution for 
one of the 77 baseline sample events can be found in 
Fig. [ini The method is validated in MC samples with 
and without intrinsic Figure [T71 demonstrates that 
the METMODEL Successfully predicts the shape of 
distributions in pythia 77 and pythia Z-^e+e^ events 
with fake frr- The technique is also cross-checked by 
performing the entire analysis with the data 77 control 
sample and data Z^e~^e~ sample. 

The V{f}r) distribution for a given event can be ob- 
tained using a large number of pseudoexperiments. For 
each pseudo-experiment, we start by forming a list of all 
jets with i?T>3 GeV and |77|<3.0 in this event and then 
smear their energies according to JER{E^^'^ ,r]) described 
above. If the smeared jet energy, is above the 

15 GeV threshold, we calculate the contribution of that 



formula: 



jet,i 



jet into the fake pr- Pr =Et - E'jP^''''\ Therefore, 
the METMODEL should account for a correlation between 
the directions of fr and jets. Then, we re-calculate the 
unclustered energy based on E^^"""^ of all jets to avoid 
double-counting when one of the jets with i?T<15 GeV 
has _E|.""^"'">15 GeV. For the next step, we randomly 
generate the expected x and y components of the frr con- 
tribution due to the unclustered energy deposited in the 
calorimeter. This procedure also accounts, on average, 
for effects of energy resolution of photons and electrons 
as well as residual effects of the wrong vertex choice. Fi- 
nally, we take a vector sum of all individual fhr compo- 
nents due to the soft unclustered energy and each of the 
jets with E^'^'^'^>lb GeV to obtain the final prediction 
of the fake fr- 

The METMODEL is not designed to predict the exact 
value of the fake fr in each event. Instead, it provides 
a two-dimensional probability density function, Vifr), 
for values of the fake fir which could arise from energy 
mismeasurements in the calorimeter. This V^fr) can be 
used to determine a significance of the observed fr in a 
given event according to the following formula: 



-significance — — log 



10 



1 - 



V{z)dz] , (CI) 



where z is the generated fake fr and w is the observed 
fjT- The .^-significance defined by Eq. ICll correctly 
takes into account all of the correlations between jets and 
the observed fir- However, the method has one signifi- 
cant drawback since it requires generating a large number 
of pseudoexperiments (e.g., >10^ pseudoexperiments for 
_^-significance=6). To overcome this problem, we take 
a simplified path of calculating an upper limit on the 
^^-significance ( "raw" .^-significance) according to the 



r-significance = -logioiVjetsVund), (C2) 

Vuncl = n f 1 - / ' KncMdu) , 



— IP IP y 



7 T^-i : 



jets 



i,Vi>0 



jets 



Vjets = n / ^'(")^" X n ^^(«)'^« 



i,t)i<0 



V, ^ f}r/{EjrCosA(l)i), 



where 'P^l^^i{u) is the probability density function for 
unclustered energy contribution to fr resolution (illus- 
trated in FigfHk). Vi{v) is the probability density func- 
tion for jet energy resolution (shown in Fig fT5|) . E'^ 
is the transverse energy of the i-th jet, and cosA(f)i 
is the azimuthal angle between that jet and measured 
fr- The "raw" .^-significance obtained from Eg. IC2l is 
then calibrated to have a simple shape defined a pri- 
ori: diV/da;=iVe«„fln(10.0)-10~^, where x is the fr- 
significance and N^vnt is the number of events in a sam- 
ple. The shape of the .^-significance has one important 
property: if all events in a data sample were to have 
only fake fr, then Nevnt ■ 10"™* events would pass a 
requirement .^-significance>cut. This property makes 
it very easy to calibrate the .^-significance by means 
of pseudoexperiments. In each pseudo-experiment, we 
obtain a randomly generated value of fr- Then we 
calculate the significance of this generated fr as if it 
were measured fr- We repeat this procedure for all 
events in the data sample and obtain the significance 
distribution for pseudoexperiments. Finally, an adjust- 
ment factor is derived for each bin of the distribu- 
tion so that the corrected .^-significance satisfies the 
A^(.^-significance>CMt)=iVei,„t • 10~™* requirement. 

The systematic uncertainties associated with the met- 
MODEL predictions are evaluated by comparing the re- 
sults obtained with the default set of parameters to 
predictions obtained with the metmodel parameters 
changed by one standard deviation (icr). In total, ten 
sources of the systematic uncertainties are considered: 1) 
difference in the unclustered energy parameterization of 
the fr resolution for 77 control and Z— >-e+e~ events; 
2) uncertainties on four parameters of the unclustered 
energy parameterization; 3) uncertainties on five param- 
eters of the JER parameterization. The correlations be- 
tween these parameters are also taken into account. The 
statistical uncertainty that depends on the number of 
pseudoexperiments per event and the systematic uncer- 
tainty are added in quadrature to obtain the total uncer- 
tainty. 
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